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SUMMARY 


This  report  describes  the  results  of  a  research  effort  to  in¬ 
vestigate  a  method  of  obtaining  high  resolution  images  of  space  ob¬ 
jects  using  earth-bound  optical  telescopes  despite  the  turbulence  of 
the  atmosphere.  The  results  of  this  research  are  an  indication  that, 
using  an  iterative  reconstruction  algorithm,  it  is  feasible  to  re¬ 
construct  diffraction-limited  images  from  the  Fourier  modulus  (or 
autocorrelation)  data  provided  by  stellar  speckle  interferometry. 

Experiments  were  performed  on  astronomical  data.  It  was  neces¬ 
sary  to  develop  methods  of  compensating  for  systematic  errors  and 
noise  in  the  data.  These  methods  were  applied  to  binary  star  data, 
and  a  diffraction-limited  image  was  successfully  reconstructed  from 
the  resulting  Fourier  modulus  data. 

The  uniqueness  of  images  reconstructed  from  Fourier  modulus  data 
was  explored  using  the  theory  of  analytic  functions.  It  was  shown, 
among  other  things,  that  if  an  object  or  its  autocorrelation  consists 
of  separated  parts  satisfying  certain  disconnection  conditions,  then 
it  is  usually  uniquely  specified  by  its  Fourier  modulus. 

A  new  method  was  developed  for  reconstructing  the  support  of  an 
object  from  the  support  of  its  autocorrelation;  it  involves  taking 
the  intersection  of  three  translates  of  the  autocorrelation  support. 
For  objects  consisting  of  a  number  of  separated  points,  a  new  method 
was  developed  for  reconstructing  the  object. 
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FOREWORD 

Tins  report  was  prepared  by  the  Radar  and  Optics  Division  of  the 
Environmental  Research  Institute  of  Michigan.  The  work  was  sponsored 
by  the  Air  Force  Office  of  Scientific  Researcn/AFSC,  United  States 
Air  Force,  under  Contract  No.  F49620-80-C-0006. 

This  interim  scientific  report  covers  work  performed  between  1 
October  1979  and  30  September  1980.  The  contract  monitor  is  Dr. 
Henry  Rauoski,  Directorate  of  Physical  and  Geophysical  Sciences, 
AFUSR/NP,  Building  410,  Bolling  Air  Force  Base,  D.C.  20332.  The 
principal  investigator  is  James  R.  Fienup.  Major  contributors  to 
the  effort  are  Thomas  R.  Crimmins  and  James  R.  Fienup.  Additional 
contributors  to  the  effort  are  Gerald  B.  Feldkamp,  Lawrence  S.  Joyce, 
Emmett  N .  Leith,  and  Christopher  J.  Roussi. 
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HIGH  RESOLUTION  IMAGING  OF  SPACE  OBJECTS 


1 

INTRODUCTION  AND  OBJECTIVES 

This  report  describes  the  results  of  the  first  year  of  a  two-year 
research  effort  to  investigate  a  method  of  ootaining  high  resolution 
images  of  space  objects  using  earth-bound  optical  telescopes. 

A  serious  problem  in  astronomy  is  that  the  turbulence  of  the 
earth's  atmosphere  severely  limits  the  resolution  of  large  earth- 
bound  optical  telescopes.  Under  good  "seeing"  conditions  the  reso¬ 
lution  allowed  by  the  atmosphere  is  typically  one  second  of  arc,  com¬ 
pared  with  0.02  seconds  of  arc,  the  theoretical  diffraction-limited 
resolution  of  a  five-meter  diameter  telescope.  That  is,  the  poten¬ 
tial  exists  for  obtaining  images  having  fifty  times  finer  resolution 
than  what  is  ordinarily  obtainable. 

Several  interferometric  methods  are  capable  of  providing  high- 
resolution  (diffraction-limited)  information  through  atmospheric  tur¬ 
bulence.  The  most  promising  of  these  interferometric  methods  is 
Labeyrie's  stellar  speckle  interferometry.  The  high-resolution  in¬ 
formation  provided  by  these  methods  is  the  modulus  of  the  Fourier 
transform  of  the  object;  the  phase  of  the  Fourier  transform  is  lost. 
Unfortunately,  except  for  the  very  special  case  in  which  an  un¬ 
resolved  star  is  very  near  the  object  of  interest,  the  Fourier  modu¬ 
lus  can  be  used  to  directly  compute  only  the  autocorrelation  of  the 
object,  but  not  the  object  itself.  The  autocorrelation  is  ordinarily 
useful  only  for  determining  the  diameter  of  the  object  or  the  sepa¬ 
ration  of  a  binary  star  pair. 

In  recent  years  it  has  been  shown  that  this  stumbling  block  can 
be  overcome  by  an  iterative  algorithm  for  computing  the  object's 
spatial  (or  angular)  brightness  distribution  from  its  Fourier  modu¬ 
lus.  The  algorithm  relies  both  on  the  Fourier  modulus  data  measured 
by  stellar  speckle  interferometry  and  on  the  a  priori  constraint  that 
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the  object  distribution  is  a  nonnegative  function.  Therefore,  the 
combination  of  stellar  speckle  interferometry  with  the  iterative 
algorithm  can  provide  diffraction-limited  images  despite  the  presence 
of  atmospheric  turbulence. 

The  goals  of  this  two-year  research  effort  are  threefold:  (1)  to 
improve  the  iterative  reconstruction  algorithm  to  make  it  operate 
reliably  in  near  real-time  on  imperfect  real-world  data,  (2)  to  de¬ 
termine  the  uniqueness  of  the  solution  under  various  conditions,  and 
(3)  to  demonstrate  the  reconstruction  technique  with  real-world 
interferometer  data,  thereby  providing  images  with  finer  resolution 
than  would  ordinarily  be  possible. 

As  envisioned  in  the  statement  of  work  for  this  contract,  these 
goals  would  be  met  as  follows: 

A.  Perform  initial  studies  and  set  priorities  for  the  following 

f i ve  study  areas : 

1.  analytical  study  of  the  input-output  concept  using  a 
statistical  approach. 

2.  analytical  and  computer  studies  of  the  uniqueness 
problem. 

3.  variations  of  the  basic  algorithm  to  improve  reliability. 

4.  analysis  and  computer  simulations  of  the  effects  of  noise 
and  other  imperfections  in  the  data,  and  methods  for 
minimizing  their  effects  for  the  types  of  noise  present 
in  conventional  interferometers. 

b .  combining  the  iterative  approach  with  other  imaging 
techniques  such  as  the  Knox-Thompson  method. 

B.  Perform  detailed  studies  of  those  areas  listed  above  that 

are  found  to  be  most  important. 
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C.  Obtain  interferometer  data,  evaluate  it,  and  process  it  into 
imagery . 

D.  Study  the  applicability  of  the  iterative  technique  to  other 
problems . 
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2 

RESEARCH  ACCOMPLISHMENTS 


Trie  first  year's  research  effort  can  be  divided  into  three  major 
topics.  (1)  Stellar  speckle  interferometer  data  was  acquired,  eval¬ 
uated,  anu  processed  into  imagery.  Methods  were  developed  for  mini¬ 
mizing  the  effects  of  the  types  of  noise  and  imperfections  found  in 
that  uata'.  (2)  Analytical  (and  to  a  lesser  extent  computer) 
studies  of  tne  uniqueness  problem  were  performed  ’  .  (3)  A  new 
method,  not  envisioned  at  the  beginning  of  this  program,  was  devel¬ 
oped  for  reconstruct ing  the  support  of  an  object;  and  for  objects 
consisting  of  a  number  of  point-like  sources,  a  new  noniterative 
method  was  developed  for  reconstructing  the  object^.  Publications 
arising  from  this  research  effort  are  listed  as  References  1-4  below. 

The  results  obtained  for  each  of  the  three  topics  listed  above 
are  summarized  in  the  three  respective  sections  that  follow.  Refer¬ 
ence  1  to  4  are  included  as  Appendices  A,  B,  C,  and  D,  respectively. 
Chapter  3  of  this  report  contains  a  more  complete  discussion  of  the 
support  reconstruction  method. 


1.  J.R.  Fienup  and  G.B.  Feldkamp,  "Astronomical  Imaging  by  Process¬ 
ing  Stellar  Speckle  Interferometry  Data,"  presented  at  the  24th 
Annual  Technical  Symposium  of  the  SPIE,  San  Diego,  Calif.,  30 
July  1981);  and  published  in  SPIE  Proceedings  Vol.  243,  App 1  ica- 
tions  of  Speckle  Phenomena  (July  1980),  p.  95. 

2.  T.R.  Crimmins  and  J.R.  Fienup,  "Phase  Retrieval  for  Functions 
with  Disconnected  Support,"  submitted  to  J.  Math  Physics. 

3.  T.R.  Crimmins  and  J.R.  Fienup,  "Comments  on  Claims  Concerning 
the  Uniqueness  of  Solutions  to  the  Phase  Retrieval  Problem," 
submitted  to  J.  Opt.  Soc.  Am. 

4.  J.R.  Fienup  and  T.R.  Crimmins,  "Determining  the  Support  of  an 
Object  from  the  Support  of  Its  Autocorrelation,"  presented  at 
the  1980  Annual  Meeting  of  the  Optical  Society  of  America, 
Chicago,  Ill.,  15  October  1980;  Abstract:  J.  Opt.  Soc.  Am.  70, 
1581  (1980). 
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<M  ASTRONOMICAL  DATA  PROCESSING 


Stellar  speckle  interferometry  data  was  obtained  both  from  the 
Steward  Observatory  Stellar  Speckle  Interferometry  Program  (via  K. 
Hege,  Steward  Observatory)  and  from  the  Anglo-Austra 1 ian  Telescope 
(via  J.C.  Dainty,  U.  Rochester). 

The  data  from  the  Anglo-Austral ian  Telescope  was  in  the  form  of 
many  short-exposure  images  on  16  mm  cine  film.  A  number  of  methods 
for  digitizing  the  data  were  explored,  and  the  one  chosen,  the  most 
economical  by  far,  was  the  following.  The  16  mm  film  was  cut  into 
strips  and  contact  copied,  along  with  grey-scale  step  wedges,  onto  9 
x  9  inch  sheets  of  film.  The  9-inch  sheets  of  film  were  then  sent 
to  the  Image  Processing  Institute  at  the  University  of  Southern 
California  for  digitization  on  their  Optronics  digitizer.  Software 
would  then  have  to  be  developed  in  order  to  extract  the  desired  data 
from  the  digitized  array  (which  includes  16  mm  film  sprocket  holes, 
etc.).  After  the  data  was  digitized  it  was  discovered  that  the 
Optronics  digitizer  had  been  malfunctioning  and  required  repairs. 
Ine  film  will  have  to  be  redigitized  before  further  experimentation 
can  proceed  with  this  data. 

Considerable  progress  was  made  with  the  Steward  Observatory  data, 
which  was  already  in  digital  form.  A  description  of  that  work  is 
found  in  Appendix  A  (Ref.  1),  and  is  summarized  below. 

It  was  previously  known  that  it  is  necessary  to  compensate  the 
Fourier  modulus  data  for  a  certain  noise  bias  term  due  to  photon 
noise.  Using  the  Steward  Observatory  data,  it  was  found  that  the 
detection  process  resulted  in  a  frequency  transfer  function,  which 
we  call  the  detection  transfer  function,  which,  in  addition  to  being 
an  error  itself,  prevented  the  compensation  of  the  noise  bias. 
Methods  of  determining  the  detection  transfer  function  from  the  data 
arid  compensating  for  it  was  developed.  Methods  of  compensating  for 
other  systematic  errors  were  also  developed.  These  methods  were 


applied  to  stellar  speckle  interferometry  data  of  a  binary  star  sys¬ 
tem,  and  a  diff ract ion-1 imi ted  imaye  was  successfully  reconstructed 
from  tin?  resulting  compensated  Fourier  modulus  data. 

Having  gained  this  experience  with  single  and  binary  star  data, 
toe  next  step  will  be  to  use  the  same  methods  on  more  complicated 
objects,  such  as  asteroids  or  Jovian  moons. 

.'.o  UNlqUtNESS  THtUKY 

Tne  principle  means  of  exploring  the  uniqueness  of  images  recon¬ 
structed  from  Fourier  modulus  data  has  been  the  theory  of  analytic 
functions.  As  described  in  more  detail  in  Appendix  B,  for  the  one¬ 
dimensional  case  there  are  usually  many  different  objects  having  the 
saii.u  Fourier  modulus.  Examples  of  both  uniqueness  and  non-uniqueness 
are  given.  However,  it  is  shewn  ttiat  if  a  function  or  its  auto¬ 
correlation  satisfy  certain  disconnection  conditions,  then  the  solu¬ 
tion  is  unique  unless  the  separated  parts  of  the  function  are  related 
to  one  another  in  a  special  way.  Therefore,  a  functions  satisfying 
these  conditions  can  usually  be  uniquely  reconstructed  from  its 
ruurier  modulus  (or  from  its  autocorrelat ion) .  It  is  also  shown  that 
it  tne  nun-real  complex  zeroes  of  the  Fourier  transform  of  a  function 
of  disconnected  support  are  finite  in  number,  then  the  support  of 
the  function  as  well  as  the  function  itself  satisfy  some  special 
.mentions.  This  makes  it  unlikely  that  the  Fourier  transform  of  a 
given  function  would  have  only  a  finite  number  of  non-real  zeroes. 

In  the  course  of  this  work  it  was  aiscovereo  that  some  of  the 
theory  appearing  previously  in  the  literature  was  in  error,  as  de¬ 
scribee  in  Appendix  C  (Ref.  3).  The  detailed  corrected  theory  is 
contained  in  Appendix  B. 

Iu  date  the  theory  of  analytic  functions  has  not  been  extended 
to  two  dimensions,  the  case  of  most  interest  in  this  research  effort. 
The  two  dimensional  case  is  not  a  direct  extension  of  the  one- 
dimensional  analysis.  The  high  probability  of  ambiguous  solutions 
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in  one  dimension  does  not  seem  to  be  the  case  in  two  dimensions.  In 
a  one-dimensional  computer  experiment  using  the  iterative 
reconstruct  ion  algorithm  on  a  case  known  to  have  two  solutions 
(Figure  1  of  Appendix  B),  the  algorithm  converged  to  one  of  the  so¬ 
lutions  in  about  half  of  the  trials  and  converged  to  the  other  solu¬ 
tion  in  the  other  half  of  the  trials,  depending  on  the  random  number 
sequence  used  as  the  initial  input  to  the  algorithm.  Therefore,  it 
is  believed  that  if  there  are  multiple  solutions,  then  the  algorithm 
is  likely  to  find  any  one  of  them.  For  the  case  of  complicated  two- 
dimensional  objects  on  the  other  hand,  the  algorithm  generally  con¬ 
verges  to  the  object  itself,  and  not  to  other  solutions.  This  is  an 
indication  that  other  solutions  do  not  usually  exist  in  the  two- 
dimensional  case. 

2.3  NEW  METHODS  FOR  SUPPORT  AND  OBJECT  RECONSTRUCTION 

As  discussed  in  more  detail  in  Chapter  3  of  this  report  (see  also 
Appendix  0,  Ref.  4),  a  new  method  was  developed  for  reconstructing 
the  support  of  an  object  (the  set  of  points  at  which  it  is  nonzero) 
from  the  support  of  its  autocorrelation.  In  some  instances,  for  ex¬ 
ample  to  find  the  relative  locations  of  a  collection  of  point-like 
stars,  the  object's  support  is  the  desired  information.  More  gener¬ 
ally,  once  the  object's  support  is  known,  then  the  complete  re¬ 
construction  of  the  object  by  the  iterative  method  is  simplified. 
Several  methods  are  shown  of  finding  sets  which  contain  all  possible 
support  solutions.  Particularly  small  and  informative  sets  contain¬ 
ing  the  solutions  are  given  by  the  intersections  of  two  translates 
of  the  autocorrelation  support.  For  the  special  case  of  convex  ob¬ 
jects,  the  intersection  of  three  translates  of  the  autocorrelation 
support  generates  a  family  of  solutions  to  the  support  of  the  object. 
For  the  special  case  of  an  object  consisting  of  a  collection  of 
points  satisfying  certain  nonredundancy  conditions,  the  intersection 
of  three  translates  of  the  autocorrelation  support  generates  a  unique 


15 


solution.  In  addition,  for  these  same  objects,  by  taking  the  product 
ot  three  translates  of  the  autocorrelation  function,  one  can  recon¬ 
struct  the  object  itself  in  addition  to  reconstructing  the  support 
ot  the  object. 

L\4  CUNILUMUNS 

All  ot  the  results  noted  above  are  encouraging  and  are  further 
indications  that  nigh  resolution  imaging  by  combining  the  iterative 
algorithm  with  stellar  interferometry  data  is  feasible.  The  pre¬ 
liminary  experience  with  astronomical  data  shows  that  although  addi¬ 
tional  problems  exist  with  real-world  data,  the  problems  encountered 
so  tar  can  be  overcome,  and  it  is  possible  to  reconstruct  high- 
resolution  hnages  from  such  data.  Fears  that  the  Fourier  modulus 
data  might  admit  to  multiple  image  solutions  are  largely  unjustified. 
The  theory  of  analytic  functions  predicts  that  a  large  class  of  ore- 
uimensiorial  functions  are  uniquely  specified  by  their  Fourier  modu¬ 
lus;  in  addition,  for  the  more  practical  two-dimensional  case  it 
appears  that  the  vast  majority  of  functions  are  uniquely  determined 
oy  the  Fourier  modulus.  Finally,  new  methods  were  oeveloped  for  re¬ 
constructing  the  object's  support  from  its  autocorrelation's  support, 
and  even  for  reconstructing  the  object  itself  by  a  very  simple  method 
for  the  case  ot  a  collection  of  point-like  stars. 
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RECONSTRUCTION  OF  THE  SUPPORT  OF  AN  OBJECT 
FROM  THE  SUPPORT  OF  ITS  AUTOCORRELATION 

3.1  I NTROOUCTIUN 

In  astronomy.  X-ray  crystallography  and  other  disciplines  one 
often  wishes  to  reconstruct  an  object  from  its  autocorrelation  or, 
equivalently,  from  the  modulus  of  its  Fourier  transform  (i.e.,  the 
phase  retrieval  problem)  .  It  is  also  useful  to  be  able  to  re¬ 
construct  just  the  support  of  the  object  (the  set  of  points  over 
which  it  is  nonzero).  In  some  cases,  for  example,  to  find  the  rela¬ 
tive  locations  of  a  number  of  point-like  stars,  the  object's  support 
is  the  desired  information.  In  addition,  once  the  object's  support 
is  known,  the  reconstruction  of  the  object  by  the  iterative 
method^  is  simplified.  Therefore,  we  are  motivated  to  find  a 
quick  way  to  determine  the  support  of  the  object  from  the  support  of 
its  autocorrelation. 

In  the  general  case  there  may  be  many  solutions  for  the  object's 
support  given  the  autocorrelation  support.  In  what  follows  a  method 
for  generating  sets  containing  all  possible  solutions  is  given.  In 
addition,  for  the  special  case  of  convex  sets  a  method  for  generating 
a  family  of  support  solutions  is  described.  For  the  special  case  of 
point-like  objects  this  method  is  shown  to  yield  a  unique  support 
solution  unless  the  vector  separations  of  the  points  in  the  object 
satisfy  certain  redundancy-type  conditions.  If  instead  of  manipu¬ 
lating  the  autocorrelat ion  support  one  uses  the  autocorrelation 
function,  then  for  the  same  point-like  objects  one  can  reconstruct 
the  object  itself.  In  the  following,  several  lengthy  proofs  are 
omitted  for  the  sake  of  brevity. 


5.  See,  for  example,  H.P.  Baltes,  ed..  Inverse  Source  Problems  in 
Optics  ( Spr inger-Ver lag.  New  York,  19787^ 

6.  J.R.  Fienup,  "Reconstruction  of  an  Object  from  the  Modulus  of 
Its  Fourier  Transform,"  Opt.  Lett.  _3»  27  (1978);  J.R.  Fienup, 
"Space  Object  Imaging  Through  the  Turbulent  Atmosphere,"  Opt. 
Eng.  J8,  529  (1979). 
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s.t  ULF  1N1  I  IllNb  AND  BACKGROUND 


Ihe  results  shown  here  apply  to  functions  of  any  number  of  di¬ 
mensions  except  where  otherwise  noted.  For  simplicity  we  consider 

only  real,  norineyative  functions.  A  function  f ( x )  >  0,  where 

N 

XfK  ,  has  support  S,  where  S  is  the  smallest  closed  set  outside  of 
which  the  function  is  ^ero  almost  everywhere. 

The  aut ucorre 1  at i on  of  f (x)  is 


t*f(x)  =  /■••/  f(y)  f(y +  x)  dV(y)  u) 

—  30 

=  /•••/  f(y)  f(y  -  x)  dv(y)  (2) 

N 

wiiei'e  V  is  the  volume  measure  on  E  .  The  Fourier  transform  of  the 
autocurre  i -L  ion  of  f(x)  is  equal  to  the  squared  modulus  of  the 
Fourier  transform  of  f(x).  Note  that  the  autocorre lat ion  is 
Ueiitro-)  symmetric.  It  is  most  illuminating  to  interpret  Eq.  (2) 
as  a  weighted  sum  of  translated  versions  of  f(x).  That  is,  in  the 
integrand  of  Eq.  (2a),  f(y)  acts  as  the  weighting  factor  for  f(y  + 
x),  which  is  f ( x )  translated  by  -y.  [t  can  be  shown  that  the  support 
of  the  autocorrelation  of  f(x)  is 


A  -  ■  '  0>  y) 

yeS 

=  b-  b-  jx-y:  x,  y  c  b  (  (3) 

Note  that  A  is  symmetr ic : 

-A  -  A  (4) 


To  illustrate  the  interpretation  of  an  autocorrelation  support,  con¬ 
sider  the  case  of  the  two-dimensional  support  S  shown  in  Figure  1(a), 
having  the  form  of  a  triangle  with  vertices  at  points  a,  b,  and  c. 
The  autocorrelation  support  A  can  be  thought  of  as  being  formed  by 
successively  translating  S  so  that  each  point  in  S  is  at  the  origin, 
and  taking  the  union  of  all  these  translates  of  S.  Figure  1(b) 
shows  three  such  translates,  (S  -  a),(S  -  b),  and  (S  -  c).  The  rest 
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of  A  is  filled  in,  as  shown  in  Figure  1(c),  by  including  all  (S  -  y) 
such  that  y  t  b. 

we  are  concerned  with  t fie  following  problem.  Given  a  symmetric 
set  A  t ^  f  i no  sets  b  C  E^  which  satisfy  A  =  S  -  S. 

bets  j  |  and  b,  are  equivalent, 

S j  ~  b^  ( ba ) 

it  there  exists  a  point  v  such  that 

b^  =  v  +  8b |  (bb) 

=  {v  +  ex  :  x  E  b1  }  (5c) 

where  a  =  +  1  or  -  l .  From  Eq.  3(b)  it  is  easily  seen  that  if  b| 
ij  a  solution  to  b  -  b  =  A,  and  if  5^  ~  b^,  then  is  also  a 
solution.  If  b|  is  a  solution  and  all  other  solutions  are  of  the 
form  v  r  ^5  then  the  solution  is  said  to  be  unique  and  A  is  said 
to  be  unambiguous;  if  there  exist  any  nonequivalent  solutions,  then 
the  solution  is  nonuri  ique  and  A  is  ambiguous . 

Not  all  symmetric  sets  are  necessarily  autocorrelation  supports. 
For  non-null  sets  A  it  follows  from  Eq.  (3)  that 

0  c  A  (6) 

is  j  necessary  condition  for  the  existence  of  a  solution.  The  fol¬ 
lowing  example  shows  that  this  is  not  a  sufficient  condition.  As 
shown  in  Figure  b,  let  A  =  )  (0,  0),  (1,  0),  (-1,  0),  (0,  1),  (0, 
-1](.  because  of  the  points  (0,  0),  (1,  0),  (-1,  0)  a  solution  must 
include  two  points  separated  by  (1,  0).  Similarly  because  of  the 
points  (0,  0),  (0,  I),  (0,  -1)  a  solution  must  include  two  points 
separated  by  (0,  1).  Therefore,  the  solution  must  have  at  least 
three  distinct  noncolinear  points.  Of  the  three  possible  pairings 
of  the  three  points,  one  has  a  separation  along  (1,  0),  a  second  has 
a  separation  along  (0,  1),  and  the  third  pair  of  points  must  have  a 
diagonal  separation.  However,  no  diagonal  terms  appear  in  A,  and 
therefore  there  is  no  solution  for  A  =  5  -  b  in  this  case. 


19 


S-b 


(b) 


(c) 


Figure  1.  Autocorrelation  Support,  (a).  Set  S;  (b)  three  of 
the  translates  of  S  that  make  up  A;  (c)  autocorrelation  support 

A  =  S-S. 


Figure  2.  A  Symmetric  Set  that  Is  Not 
an  Autocorrelation  Support 
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A  set  X  is  convex  if  for  all  x,  y  e  X, 

tx  Ml  -  t)y  e  X  (7) 

f or  al 1  t  e  [0,  1 ] .  The  convex  hu 1 1  of  a  set  X,  denoted  by  c.hul 1 (X ) 

N 

is  given  by  the  smallest  convex  subset  of  E  containing  X.  Thus  X 

is  convex  if  and  only  if  X  =  c.hull(X).  If  S  is  convex,  then  A  =  S 

-  3  is  also  convex.  More  generally, 

c.hul 1(X  -  X)  =  c.hul 1(X)  —  c.hull(X)  (8) 

All  convex  symmetric  sets  A  have  at  least  one  solution 

S  =  1  A  =  {  x/2  :  x  e  A  |  (9) 

Trie  proof  is  as  follows.  Let  u,  v  e  1/2  A.  Then  2u  e  A,  2v  e  A  and 
-2 v  r  A.  Therefore, 

u  -  v  =  ^  (2u )  ♦  \  {-2m)  e  A  (10) 

arid  so  (1/2  A)  -  (1/2  A)  C  a.  NOW  let  v  E  A.  Then  v/2  e  1/2  A  and 

-v/2  e  1/2  A.  Therefore 

v  =  (^v)  -  (-  }v)  e  (]a)  -  (^A)  (11) 

arid  so  A  Er  (1/2  A)  -  (1/2  A).  Therefore 

A  =  (±A)  -  (^A)  (12 


3.3  LOCATOR  SETS 

In  many  cases  A  is  ambiguous,  and  it  would  be  useful  to  define  a 

N 

set  that  contains  all  possible  solutions.  A  set  L  Sr  E  is  de¬ 
fined  as  a  locator  set  for  A  if  for  every  closed  set  S  9  E^ 
satisfying  A  =  S  -  S,  some  translate  of  S  is  a  subset  of  L,  i.e., 
there  exists  a  vector  v  such  that 

v  ♦  S  9  L  (13) 
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1 hi  i  r  ,ii  i'  man y  ways  to  generate  locator  sets.  For  example,  for  v  e 
,  :  v  S  -  6  -  A ,  and  so  A  itself  is  a  locator  set.  Naturally, 

tin-  s.iM  1 lei  the  locator  set,  t he  more  tightly  it  bounds  the  possible 
s i > !  i i  ions,  ano  the  more  informative  it  is.  It  can  be  shown  that  a 
i.i  H  r  locator  set  than  A  is 

L  =  A  ■ ■  H  (14) 

.v  n  t  r  ■  t  o  is  a  closed  naif -space  with  the  origin  on  its  boundary.  A 
•ili  smaller  locator  set  can  be  shown  to  be 

L  =  l  P  (lb) 

•  e  is  any  N  -d imens iona  1  paral  lelopiped  containing  A. 

particularly  interesting  locator  set  is  given  by  the  following 
n  tersei.  f  ion  of  two  autocorre  1  at  ion  supports.  If  w  e  A,  then 

L  =  A  (w  +  A)  (16) 

is  i  locator  set  for  A.  Note  that  L  is  symmetric  about  the  point 

v  .  Ine  proot  that  this  is  a  locator  set  is  as  follows.  Suppose  S 

i  i,t  o  '  -  o  -  S.  Since  w  t  A,  there  exist  u,  v  e  S  such  that  w 

,  ..  un.iier  /  ,  S  -  v.  Then  z  =  s  -  v  whei  :  s  e  S,  z  =  s  -  v 

,  in  j  /  ■  u  <•  (u  ■  v )  ^  s  -  u  +  w  e  A  +  w.  Therefore,  z  e  A 

♦  w ;  l  and  S  -  V  t  1.  . 

Nat  .rally,  the  most  interesting  (smallest)  locator  sets  generated 
■  ,  ini,  met.  nmi  at  intersecting  two  autocorre  1  a  t  i  on  supports  are  ob- 

t.  1 1 ne. j  iiy  ,  housing  w  to  be  on  t fie  boundary  of  A.  By  choosing  dif- 

i *  •*  en  t  points  w  i  A,  a  whole  family  of  locator  sets  can  be  generated 
ay  this  HiettiOU. 

i  <  imp  le  1  . 

((insider  the  set  S  shown  in  Figure  3(a),  consisting  of  two  balls 
p  1 1  rn-ii  by  two  thin  rods,  arid  its  autocorre  1  a  t  ion  support  A  =  S  -  S 
shown  in  figure  3(b).  An  example  of  a  locator  set  1/3  P  is  shown  in 
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(f) 


Figure  3.  Locator  Sets,  (a)  Set  5;  (b)  A  =  S-S;  (c)  locator  set 
L  =  1/2  P;  (d)  formation  of  L  =  AH(w  +  A);  (e)  and  (f)  two  other 
members  of  the  family  of  locator  sets. 
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1  iqure  .He);  it  Joes  a  good  job  of  defining  the  approximate  size  of 
j,  biit  it  is  not  suggestive  of  any  of  the  details  of  its  shape. 
I  i guru  J(J)  shows  the  generation  of  the  locator  set  L  =  A  n  (w  +  A) 
tor  a  particular  point  w  c  A.  Figure  3(e)  and  3(f)  show  two  other 
memuers  of  the  family  of  locator  sets  generated  with  two  other  points 
w  t  A.  These  locator  sets  gene'ated  by  intersecting  two  auto¬ 
correlation  supports  are  very  suggestive  of  the  shape  of  the  solution 
(or-  solutions).  This  is  especially  true  if  one  realizes  that  any 
solution  must  be  contained  within  all  of  these  locator  sets.  Unfor¬ 
tunately,  for  the  general  case  it  is  difficult  to  narrow  down  the 
solution  any  further;  a  way  to  combine  the  information  from  two  or 
more  of  the  family  of  locator  sets  has  not  been  devised.  However, 
as  will  oe  shown  in  the  sections  that  follow,  for  special  classes  of 
sets  much  more  can  be  done. 

f.4  AuTJCUkRELATION  TRI-INTERSECTION  FOR  CONVEX  SETS 

Fur  the  special  case  of  convex  sets  A,  a  family  of  solutions  can 
oe  generated  by  a  simple  method.  For  the  one-dimensional  case  the 
result  is  trivial:  a  unique  solution  is  given  by  S  =  1/2  A,  which  is 
just  a  segment  of  the  line  half  the  length  of  the  l'ine  segment  A. 
An  equivalent  result  for  tne  one-dimensional  convex  case  is  the 
Su  I ut.  ion 

s  =  A  n  (w  ♦  A)  (17) 

where  w  is  on  the  boundary  of  A  (at  one  end  of  the  line  segment  A), 
or  in  symbols  w  e  3(A). 

For  tne  two-dimensional  convex  case,  we  have  the  following  re¬ 
sult.  let  A  :  t  be  a  closed  convex  symmetric  set  with  non-null 
interior,  and  let 

w j  c  3(A)  and  w?  e  3(A)  f  3(w-|  +  A). 
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Furthermore,  let 

B  =  A  i  ;  ( w -j  +  A)  n  (w 2  +  A).  (19) 

Then 

A  =  B  -  B.  (20) 

The  lengthy  proof  of  this  result  is  omitted  for  the  sake  of  brevity. 

tx ample  2 

Consider  the  set  S  shown  in  Figure  4(a),  which  is  the  convex  hull 
of  the  set  shown  in  Figure  3(a).  Its  autocorrelation  support  A  =  S 
-  S  (.which  is  the  convex  hull  of  Figure  3(b)]  is  shown  in  Figure 
4(b).  The  parallelogram  shown  in  Figure  3(c)  is  a  locator  set  for 
A.  A  member  of  the  family  of  locator  sets  A  H  (w  +  A)  is  shown  by 
the  intersection  of  A  and  w  +  A  in  Figure  4(c).  A  member  of  the 
family  of  solutions  B  is  shown  by  the  intersection  of  the  three  sets 
A  n  (w^  +  A)  n  (w.p  +  A)  in  Figure  4(d).  Two  other  examples  of 
b  are  shown  in  Figures  4(e)  and  4(f). 

3.6  THREE-DIMENSIONAL  INTERSECTIONS  OF  CONVEX  SETS 

For  convex  sets,  since  in  one  dimension  the  intersection  of  two 
sets,  Eq.  (17),  results  in  the  solution,  and  since  in  two  dimensions 
the  intersection  of  three  sets,  Eq.  (19).  results  in  solutions,  one 
might  hope  that  in  three  dimensions  the  set 


C  =  A  O  (w]  + 

A)  n  (w2 

+  A) 

n  (w3  +  a) 

(21) 

would  be  a 

solution  to  S  ■ 

-  S  =  A, 

where 

W1  e  3(A), 

w2  e  3(A)  n 

a(W]  +  A), 

and  w^  e  a(A) 

n  3  ( w  i 

+  A) 

n  3(w2  + 

A).  Unfortu- 

nately,  this  is  generally  not  the  case. 

A  counter-example  to  C  -  C  =  A  is  the  following.  Consider  S 
equal  to  a  sphere  of  diameter  one,  then  A  =  5  -  S  is  a  sphere  of 
radius  one  centered  at  the  origin.  Figures  5(a)  and  5(b)  show 
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Figure  4.  Autocorrelation  Tri  -  Intersection  Solutions  for  Convex  Sets 
(a)  Set  S;  (b)  A  =  S-S;  (c)  formation  of  locator  set  L  =  AO(w  +  A); 
(d)  formation  of  solution  B  =  AO(w,  +  A)n(wp  +  A);  (e)  and  (f)  two 
other  solutions  Af  the  form  B. 


Figure  5.  Sphere/Circle  Exam 
(c)  B  =  APl(w1  +  A)n(w 2  +  A) 


,)  1  an  a*'  cuts  through  the  centers  of  S  and  A,  respectively.  Figure 
->(<;)  shows  a  planar  cut  through  A  n  (wj  +  A)  r.  (w^  +  A) 

t'uough  the  tliree  points,  0,  W|  and  w^.  Art  (w^  +  A)  f\ 

•w„  +  A)  has  two  vertices,  one  in  front  of  the  plane  of  the  page 
and  one  Dehind  the  plane  of  the  page,  both  at  distance  one  from  the 
.(.•titers  of  each  of  the  three  intersecting  spheres.  Taking  the  inter¬ 
section  of  this  with  (w^  +  A),  which  is  centered  at  one  of  the  two 
vertices,  gives  us  C,  which  is  similar  to  a  regular  tetrahedron  (it 
has  toe  same  vertices)  but  having  spherical  surfaces  of  radius  one 
in  place  of  the  four  plane  faces  of  a  tetrahedron.  Looking  for  a 
inn;  i- n  t  at  the  tetrahedron  T  having  the  same  vertices  as  C  (i.e.,  the 
i.iin vc- a  hull  of  points  0,  w w^,  and  w^  having  edges  of  length 
one;,  we  see  tnat  F  -  T  is  a  cuboctahedron ,  which  has  eight  tri¬ 

angular  faces  and  six  square  faces.  Since  TC  C,  T-TCC-C. 
Tne  surface  of  C  -  C  can  be  subdivided  into  twelve  patches  associated 
with  the  twelve  faces  of  the  cuboctahedron.  It  can  be  shown  that 
the  eight  patches  associated  with  the  triangular  faces  coincide  ex- 
actly  with  the  surface  of  the  sphere  A  of  radius  one.  However,  the 
jU  patches  correspond i ng  to  the  square  faces  do  not.  For  example, 
i.ne  d  i  j  Lance  from  the  origin  to  the  center  of  each  of  those  six 

pot. i  ries  is  equal  to  the  distance  between  the  centers  of  two  opposing 
edge ,  of  i.  This  distance  can  be  shown  to  be  yjl  -yJl/2  =  1  .0249; 
L n  if  is,  L  -  C  bulges  beyond  the  sphere  by  about  2.49  percent  at 

those  points.  Therefore,  C  -  C  ^  A. 

i.b  LUMP  1  MAT  I  ON  b  UF  CONVEX  SOLUTIONS 

Solutions  to  convex  A  =  S  -  S  of  the  form  B  make  up  a  family  of 
solutions  generally  having  an  uncountable  infinity  of  members,  one 
for  each  W|  e  a(A).  Nevertheless,  there  may  exist  additional 

so  1  a t i ons  . 
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Additional  solutions  can  be  generated  in  the  following  way.  If 
b|  and  5.,  are  solutions  to  convex  A  =  S  -  S,  then 

St  =  tS,  +  (1  -  t)S2  (22) 

is  also  a  solution  for  0  <  t  <  1.  The  proof  of  this  result  is  as 

f ol lows 

St  -  St  =  Its,  +  (1  -  t)S2]  -  Its,  +  (1  -  t)S2] 

=  ts,  -  t s i  +  (1  -  t)S^  -  (1  -  t)S„ 

(23) 

=  tA  +  (1  -  t)A 
=  A 


since  A  is  convex. 


If  S |  is  a  solution,  then  so  is  -  S,.  Then  using 
ana  =  -  S,  in  Eq.  (22),  it  is  seen  that 

c  _  1  c  1  c  _  J_  a 

i 12  "  2  H  2  al  "  2  ' 


-  1/2 

(24) 


is  a  solution,  as  was  previously  shown  by  Eq.  (12). 

Eq.  (22)  can  easily  be  generalized  as  follows.  If  S , ,  ..., 

Sn  are  solutions  for  convex  A,  and  if  t , ,  ...,  tfl  >  0  and  t, 

+  t.,  +  . . .  +  tn  =  1 ,  then 

S  =  2  t.  S-  (25) 

i  =  l  1  1 


is  also  a  solution. 


Example  3 

Consider  the  two-dimensional  convex  set  S  shown  in  Figure  5(a), 
consisting  of  a  circle  of  diameter  one.  A  =  S  -  S,  consisting  of  a 
circle  of  radius  one  is  shown  in  Figure  5(d),  and  a  tri-intersection 
solution  B  is  shown  as  the  intersection  of  three  circles  in  Figure 
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5(c).  This  solution  is  similar  to  an  equilateral  triangle  but  having 
arcs  of  circles  of  radius  one  with  centers  at  the  opposite  vertices 
for  each  of  the  three  sides.  It  can  easily  be  seen  that  all  other 
solutions  b  generated  by  Eq.  (19)  are  similar  to  the  one  shown  in 
Figure  b(c)  except  rotated  in  the  plane.  The  circle  of  diameter  one 
shown  in  Figure  5(a)  is  not  of  this  form,  but  it  is  also  a  solution 
to  A.  As  shown  by  Eq.  (24),  S  =  1/2  A  in  Figure  5(a)  can  be  gene¬ 
rated  by  applying  Eq.  (22),  using  Sj  =  -  =  B  and  t  =  1/2. 

one  of  a  family  of  additional  solutions  generated  by  Eq.  (22)  is 
shown  in  Figure  5(d).  It  was  generated  using  =  1/2  A  in  Figure 
5(a),  =  B  in  Figure  5(c),  and  t  =  1/2. 

3./  THE  AMBIGUITY  OF  CONVEX  SETS 

we  now  consider  the  question  of  uniqueness  of  convex  solutions 
of  A  =  S  -  S  for  convex  A.  As  mentioned  earlier,  1/2  A  is  a  solu¬ 
tion.  If  all  convex  solutions  are  equivalent  to  S,  then  A  is  said 
to  be  convex-unambiguous .  It  was  shown  that  in  two  dimensions  one 
can  generate  a  family  of  solutions  by  Eq.  (IS),  the  member  of  the 
family  being  determined  by  the  choice  of  w^ .  Eq.  (22)  or  (25)  can 

then  be  used  to  generate  still  more  solutions.  Therefore  one  would 
suppose  that  convex  sets  A  are  generally  convex-ambiguous.  However, 
it  is  also  possible  that  all  solutions  generated  by  Eq.  (19)  are 
equivalent,  in  which  case  A  would  be  convex-unambiguous. 

In  what  follows  it  is  shown  that  in  two  dimensions  if  A  is  a 

parallelogram  then  A  is  convex-unambiguous.  Let  A  be  a  parallelo¬ 
gram  having  vertices  w^,  -  w^ ,  w^,  and  -w^ .  By  Eq.  (16)  a 

locator  set  for  A  is  L  =  A  n  (w-|  +  A)  since  w-|  e  A.  It  is  eas¬ 

ily  seen  that  L  =  1/2  w1  +  ]/2  A,  and  so  L'  =  1/2  A,  which  has 
vertices  1/2  Wp  -  1/2  w^,  1/2  w^,  -  1/2  w^  is  a  locator  set 

fur  A.  Suppose  A  =  S  -  S  where  S  is  convex.  Then  some  translate  of 

S ,  call  it  S',  is  contained  in  L'.  Since  w^  e  A  there  exist  u,  v 
e  S'  such  that  w^  =  u  -  v.  Since  S’  9  L’>  u,  v  e  L'.  It  follows 
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that  u  =  1 /2  and  v  =  -  1/2  w-j.  Therefore,  1/2  E  S'  and 
-  1/2  w |  e  S'.  Similarly  1/2  w^  e  S'  and  -  1/2  e  S'. 

Then,  since  S'  is  convex 

L'  =  c.hull  w  ^ w^,  w^,  -  j  w^jj  -  -  L‘  (26) 

Therefore,  S  =  L '  =  1/2  A,  and  so  S  is  unique  among  convex  solutions. 

It  can  also  be  shown  that  parallelograms  are  the  only  two- 
dimensional  convex-unambiguous  sets.  Convex  symmetric  sets  A  c 
V'  that  are  not  parallelograms  can  be  shown  to  have  infinitely 
many  nonequivalent  solutions  to  A  =  S  -  S. 

3.8  AUTOCORRELATION  TRI-INTERSECTION  FOR  POINT-LIKE  SETS 

For  the  special  case  of  certain  point-like  sets  A,  the  solution 
can  be  generated  by  a  method  similar  to  the  one  for  convex  sets.  By 
point-like  sets  we  mean  sets  comprised  of  a  collection  of  distinct 
noncontiguous  points.  For  example,  a  point-like  function 

N 

f(x)  =  fn  6^x  ~  xn)  (2?) 

n=l 

consisting  of  N  delta  functions  having  amplitudes  f  >  0,  n  =  1, 
...,  N,  would  have  point-like  support 

S  =  {xn  :  n  =  1,  ...,  n}.  (28) 

The  following  result  holds  for  any  number  of  dimensions.  Let  S  be  a 
point-like  set  and  A  =  S  -  S.  Let  w-j  e  A  and  w^  e  A  O  (w^  + 
A),  with  0  /  w |  ^  w^  ^  0,  and  let 

B  =  A  n  (w]  +  A)  O  (w2  +  A)  (29) 

Define  the  following  Condition  1: 
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If  V  X2 ’  yr  y2’  zr  z2  6  s>  xi  *  V  dnd 

x!  -  x2  +  y,  -  y2  +  zj  -  z2  =  0  (30) 

then 

x |  =  y^  or  z,,,  and  x2  =  or  z] 

It  can  be  shown  that  if  3  satisfies  Condition  1,  then 

B  -  S  (31) 

That  is,  B  is  the  unique  solution  to  A  =  S  -  S. 

Another  approach  is  as  follows.  Define  Condition  2:  if  the  set 
GCA  consists  of  three  distinct  points  and  if  0  e  G  and  G  -  G  C  A, 
then  G  is  equivalent  to  a  subset  of  S.  Define  Condition  3: 

if  x ( ,  x?,  y?  c  S,  Xj  +  x2,  and  Xj  -  =  y]  -  y2;  then  x]  =  y] . 

we  nave  trie  following  two  results.  If  S  satisfies  Condition  2,  then 
3  is  equivalent  to  a  subset  of  B.  It  can  also  be  shown  that  if  S 
satisfies  Conditions  2  and  3,  then  S  is  equivalent  to  B;  and  S  satis¬ 
fies  Conditions  2  and  3  if  and  only  if  it  satisfies  Condition  1. 

Since  it  requires  a  special  relationship  between  the  points  in  S 
in  order  that  Condition  1  not  be  satisfied,  it  is  probable  that  for 
3  comprised  of  -andomly  located  points,  B  is  the  unique  solution  to 
A  =  3  -  3.  More  will  ue  said  about  this  later. 

Ex ample  4 

Consider  the  point-like  set  3  having  9  points  shown  in  Figure 
b(a).  A  =  3  -  S  shown  in  Figure  6(b)  has  9^  -  9  +  1  =  73  points. 
Intersecting  A  with  a  translate  of  itself  using  Eq.  (16),  a  number 
of  different  locatur  sets  for  A  can  be  formed,  two  of  which  are  shown 
in  Figures  6(c)  and  (d).  (Any  solution  to  A  =  S  -  S  must  have  trans¬ 
lates  that  are  subsets  of  all  the  locator  sets.)  For  this  example, 
for  all  allowable  values  of  w^  and  w^,  B  is  found  to  be  equiva¬ 
lent  to  5,  which  is  shown  in  Figure  6(a). 
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It  can  also  be  shown  that  even  when  Condition  1  is  not  satisfied 
it  is  sometimes  possible  to  find  solutions  (and  the  solutions  may 
even  be  unique  as  it  was  in  Example  b)  by  intersecting  A  with  itself 
three  or  more  times.  However,  when  Condition  1  is  not  satisfied, 
then  there  is  no  guarantee  that  the  solution  is  unique,  and  finding 
solutions  is  considerably  more  complicated  than  simply  evaluating  B 
by  Lq.  ( 29 ) .  Unfortunately,  given  A  it  is  not  possible  to  immedi¬ 
ately  determine  whether  Condition  1  is  satisfied.  A  necessary  con¬ 
dition  that  Condition  1  (or  Condition  3)  be  satisfied  is  that  the 

2 

number  of  points  in  A  can  be  expressed  as  N  -  N  +  1  where  N  >  1 
is  an  integer. 

3.9  KtCUNSTKUCTION  OF  POINT-LIKE  OBJECTS 

By  a  simple  modification  of  the  method  described  in  the  previous 
section  for  reconstruct ing  the  support  of  a  point-like  object,  it  is 
often  possible  to  reconstruct  the  object  itself.  The  method  is 
analogous  to  using  Eq.  (29)  to  compute  B,  except  that  it  deals  with 
products  of  autocorrelat ion  functions  instead  of  intersections  of 
aut ocorre lat ion  supports. 

Suppose  that  the  object  is  given  by  Eq.  (27),  consisting  of  N 
delta  functions  located  at  the  distinct  points  xn  having  ampli¬ 
tudes  f  ,  n  =  1,  2,  ...,  N.  The  positions  xp  are  vectors  in  any 
number  of  dimensions.  The  autocorrelation  is 

f*f(x)  =J  f (y)  f (y  +  x)  dV(y) 

N  N 

=2  2  f  f  6(x  -  x  +  x  )  (32) 

n=l  m=l 

which  can  be  expressed  as 

/  N  \  N 

f*f(x)  =  (  X  fn)6(x)  +  X  X  ym«<*  ~  *m  +  *n)  <33) 

V  n= 1  /  n=l  m^n 


34 


\ 


which  has  terms  located  at  positions  x  =  xm  -  x  ,  N  of 

which  are  at  x  =  0.  That  is,  it  has  up  to  N?  -  N  ♦  I  distinct 
terms . 

At  this  point  we  would  like  to  take  the  product  of  two  such  auto¬ 
correlation  functions;  however,  the  product  of  two  delta  functions 
is  not  well  defined.  In  order  to  overcome  this  problem,  we  define 
tiie  product  of  two  delta  functions  as  follows: 

ab«(x  -  ) ,  x^  =  Xj 

[_a5(x  -  x  j )  ]|_bs( x  -x?)]  " 


0 


»  ^  x i 


Now  consider  multiplying  f-frf(x)  by  f*f(x  -  x,  +  xk),  where 

x,  -  x.  iu  lies  within  the  support  of  f'A'f(x).  The  center  of 
Ik 

the  translated  autocorrelation  lies  within  the  support  of  the  un¬ 
translated  autocorrelation.  This  gives  the  autocorrelat ion  product 
(all  summations  are  from  1  to  N  unless  otherwise  noted) 

AP  j R ( x )  =  If  ★  f(x)]Lf ★  f(x  -  x1  +  xk)] 

(lfn)6(x)  +  1’  X  fnfm«(x  -  xm  +  xn) 
l  n  n  m^n 


lfn)6(x  -  X1  +  V 

n 

X'  S' 

f  ,f  , 6 ( X  -  X  ,  +  X  ,  -  X,  +  X.  ) 

- n  *  in  *  111  n  1  k 


n 1  m 1  ^n 1 


n  in 


>1 


(34a) 


(34b) 


(Xfn2)flfk«(x)  +  (£  fn)flfk6^x  -  xl+  xk 


+  f  1  ^  k  4r*  .  m 
m^k ,  1 


r  «(x  -  x  -  x,  ) 

m  '  in  v  ' 


+  f,f,  X  fZ«  (x  - 
1  k  ,  n 
n^k ,  1 


x1  +  xn)  +  (O.T.; 
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rth-ro  (O.I.)  denotes  "other  terms"  as  will  be  described  later.  [As 
.in  example  of  how  tq.  34(b)  follows  from  tq.  34(a),  the  fourth  term 
i>;  iq.  qpij  arises  from  the  product  of  the  second  term  of  the  first 
.cit  in  orrv  I  at  ion  with  the  second  term  of  the  second  autocorre  I  at  ion , 
,vitn  m  =  I,  n1  =  n,  and  in'  =  k.J  From  another  way  of  expressing  tq. 
i .all 


w  N  N  N 

\  V  N  v  f  f  f  ,f  ,6(x 

. , ,  , ,  i — >  n  m  n  m 

n=  1  iii=  I  n  ’  =  1  m *  =  I 


x  +  x  ) 
m  n 


.  6  ( x  -  x  ,  +  x  ,  -  x ,  +  x  )  (34c ) 

'  m  n  1  k' 


it  i,  seen  tnat  terms  survive  at  points 


I ne  terms  shown  in  tq.  34(b)  all  necessarily  appear.  In  addition, 
,.tncr  terms  may  appear,  as  indicated  by  "+  (U.T.)".  The  existence 
I. fuer  terms  depends  on  the  presence  of  special  relationships  between 
tne  coordinates  x  allowing  Eq.  (3b)  to  be  satisfied.  There  beiny 
no  additional  terms  is  equivalent  to  Condition  1  (described  in  the 
previous  section)  being  satisfied.  If  the  x  were  independent 
rarniom  variables,  then  the  chance  of  having  additional  surviving 
terms  would  be  small,  and  we  would  have  (U.T.)  =  0. 

combining  Eq.  (27)  with  34(b),  the  autocorrelation  product  can 
be  expressed  as 

APik(x)  =  f  ifk  [f2(x  +  V  +  f2(-x  +  X1}] 


f  ,f kL & ( X )  +  5 ( x  -  x1  +  xR) j  +  (U.T.) 


(36) 


Tnerefore,  translates  of  the  supports  of  both  f(x)  and  f(-x)  are 
contained  within  the  support  of  AP^(x).  This  can  also  be  seen 
from  the  fact  that  by  Eq.  (16)  the  support  of  AP^  (x)  is  a  locator 
set . 
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Included  m  the  support  of  AP^(x)  are  points  x^  -  x^  and 

x(  x  ,  n  --  I,  Z,  ....  N.  Therefore  the  center  of  f ★  f ( x  - 
x  I  ,  +  *k),  l1  ^  1,  k  is  within  the  support  of  AP^(x).  If 

(.u.I.)  _  u,  tnen  the  product  of  the  three  autocorrelations  is 


Lt*t  (x)  J .  L  f  ★  f  (x  -  Xj  +  xk].[f*f(x  -  x  i  ,  +  xk)  ] 


,  AP|k(x)  L ★'fr ( x  -  x,  ,  +  x  J] 


1  k ; 


n  k  ' 


n^k ,  I 


Zf^5(x  -  x,  +  x  ) 
n  1  n 

n^k ,  1 


fn  6(x  -  X,.  +  xk) 


+E  E 


n'  rn'^n 


,  .  .  <5(x-x.+x,-x,,+x1) 
,  n '  m'  '  in’  n 1  1  1  k ' 


Wl 


Zf  ^6  (  X  -  X  +  X, 

n  '  n  k 

"A,  1,1* 


Ak;- 

(37a) 


,«(x) 


+  (E  f^fi6(x  -  xi 


,6(x  -  x,,  ♦  xk; 

(37b)J 


ihat  is,  the  support  of  tne  product  of  three  autocorrelations  has 
tne  s dine  support  as  f(x  +  x^),  as  was  shown  earlier  in  connection 
with  Lq.  (79),  since  B  is  just  the  support  of  the  product  of  three 
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such  autocorrelation  functions.  Furthermore,  except  at  three  points 
tne  product  is  proportional  to  the  cube  of  t(x  +  x^). 

The  values  at  all  points  can  be  determined  as  follows:  First, 

IV  =  f*f(0)  D  (38; 

n 


is  viiown,  so  that  factor  can  be  divided  out  from  the  last  three  terms 
of  Id-  3/ta).  becond,  let  the  coefficients  of  those  three  terms  in 
Ltj.  H  (a)  be  (with  X  ff  divided  out) 


n 

A  =  D" 


flpHcl-k<°>  '  <1  fl  fl 


S  .  D-1  AP,kl.K(«,  -  xk)  .  fk  f]  f, 


l  =  0  1  AP 


. u(xi.  -  =  U  fi  * 


1  k  1 '  k v  1 


( 3y  a ) 
(39b) 
(39c) 


Ivmy,  we  yet 


an  i 


1/4 


fk  = 


(40a) 


1/4 


f  _  iL. 

’  "  Vac / 


(40b) 


1/4 


f  1  ■  =  Vaiv 


(40c) 


f^y,,  =  (ABC) 


1/4 


(40d) 


1  he  remaining  t  's,  for  n  ^  k,  1,  1'  can  then  be  computed  by  di¬ 
viding  ty.  3/(a)  by  f  fjf],  and  then  taking  the  cube  root: 


f 

ri 


AP , .  ,  ,  ( x  -x  ) 
1  k  1  x  n  k 


WT 


1/3 


(40e) 
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By  In i s  method  f(x)  is  reconstructed  exactly  to  within  a  transl at  tori , 
as  long  as  (O.T. )  =  0. 

In  pertornuny  these  ca  leu  1  at  ions ,  had  we  chosen  a  translation  of 
tne  form  ( x (  -  ,  ) ,  k '  ^  k,l  instead  of  (x^,  -  ) ,  then 

the  result  would  have  been  similar,  except  a  translate  of  f(-x)  would 
have  ueen  reconstructed  instead  of  a  translate  of  f(x).  If  (O.T.)  ± 
U,  that  is,  if  Condition  1  is  not  satisfied,  then  additional  terms 
appear  that  make  the  analysis  much  more  complicated  and  may  prevent 
the  reconstruction  of  t(x). 

Various  modifications  to  this  reconstruct  ion  method  are  possible. 

For  example,  the  product  of  two  autocorrelat ion  products  AP,,  (x) 

4  K 

AP  |  , N ( * )  is  proportional  to  f  (x  +  x^)  except  at  three 
points.  Another  example  is  to  define  the  autocorrelation  support 
function  as 


N 

A  ( x )  =  6  ( x )  +V  7  6  ( x  -  x  +  x  ) 

n=l  nijtn 


(41) 


which  is  just  a  binary-valued  version  of  Eq.  (33).  Then  the  product 
uf  the  autocorrelation  function  with  two  properly  translated  auto¬ 
correlation  support  functions  is  proportional  to  a  translate  of  f(x), 
except  at.  a  single  point  which  can  be  determined  by  a  few  extra  sim¬ 
ple  steps. 


In  arriving  at  Eq.  (3/),  it  was  assumed  that  the  other  terms 
(O.T.)  -  U,  or  equivalently  that  Condition  1  be  satisfied.  The  terms 
included  in  Eqs.  34(b)  and  (37)  are  those  that  necessarily  arise  by 
satisfying 


-  x  = 

n 


+  x,  -  x, 


(42) 


trivially,  for  example,  tor  m  =  n,  m'  =  k,  and  n  ’  =  1  .  The  other 
terms  are  those  that  satisfy  Eq.  (42)  by  chance,  that  is,  those  that 
arise  in  audition  to  those  tnat  (trivially)  arise  necessarily.  These 
other  terms  require  a  special  relationship  between  the  points  in  S, 
and  would  not  be  expected  to  occur  if  the  points  in  S  are  randomly 
distributed  in  some  region  of  E^1. 
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These  results,  with  some  modifications,  can  also  be  extended  to 
the  case  of  an  object  having  support  on  a  number  of  disjoint  islands 
having  diameters  small  compared  with  their  separations  (as  opposed 
to  the  support  consisting  of  isolated  mathematical  points).  How¬ 
ever,  as  the  number  of  islands  increases  and  as  the  ratio  of  the 
diameters  ot  the  islands  to  their  separations  increases,  the  prob¬ 
ability  of  satisfying  a  condition  analogous  to  Condition  1  decreases. 
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■  :  w  *  v  H  d  >>*■  \  the-  telescope  aperture  (that  M1K  is  t.tie  autorer- 

*  f'u  t  *•  1 1'*  i  np*1  p  ip  i  1  fun  A-,suminq  a  circular  ap-rture,  a  cut  throuqh 

•  '  *•  .  r"  ■  .[  •  ip  rf'.r*  m  T  t  would  h.jv*  .j  rouqf.l/  cone  shape**  arid  tie  zero  tieynnrl  pixel  37. 

f  wi  r  ,  is*.  ..jiruru,:  pnw^r  spe-  from  <j  t  t  tie  i  j  o  i'  e  s  o  L  v  “  d  star  shown  in  figure  2  is  very  far 


:  a  .  f  d*  rr  i  i »  ,  * .  *  *  ■  e  ■  '.ape  «  f  i  ''f  pC'wer  speotrum.  First,  the  sp'dFJe  MTf^,  men- 

*.  i  i  * in  nine  t  1  on  with  fq.  (M,  drops  v-.-ry  rapidly  for  the  very  low  spatial  fre- 
po  •  s  ri'.if  re.  .j  i  'h-  r  ir.  rut -nt  r.  Itiis  result.-,  in  the  spike -like  behavior  uf  the 

-.mmi-d  power  p**»  1. 1  om  t  -i  r  very  low  spatial  frequencies.  Beyond  the  very  low  spatial- 
‘  rk  li'/'oi  ,  re;- on,  t.n»  peiKie  MTV'  is  niort*  t-etti-r  behoved  and  decreases  slowly*.  Second, 
phi. tun  e*,i  re-,  jit'  in,  amonq  other  things,  a  noise  bias  term  in  the  summed  power  spec- 
f.  r  jm  *  t  n  j  - ,  uni  •  :  ia--  term  •  ji  .m  i  na  t «  •.  in  the  fiiqtier  spatial  f  requrnr  i»*s  .  Beyond  a 

r-jdius  • f  5/  pixels  in  the  summed  power  spectrum,  no  signal  energy  exists  --  it  is  purely 
si.,  i  ,»■  .  -e  .jnm  i  n;»  *  fhe  summed  pi  wer  spectrum  that  little  useful  information  can  he 

ust  ji>»d  1 1 1 1 1  i  ■  .-mp  en  sa  t  i  on  is  made1  for  both  of  these  two  effects. 

In.  es  i  s  *•’  {.;  i  a  ,  farm  .inj  the  detei  t  ion  tl  jnsN-r  funrt  ion 


f jr --e  would  ordinarily  eliminate  the  noise  bias  term  simply  by  subtracting  a  constant  from 
t»  •  uniTed  pnw  r  ->p»  1  t  i  um  * fJ  *  ^  .  HowevfM  ,  a*.,  seen  from  Figure  7,  the  noise  1  ias  term, 

wr  is  se^Ti  fjy  itself  f*eyurid  pixel  57,  is  not.  constant  in  this  case.  This  results  from 

tr-  fet  M.jf  .rein  i  j*  - 1  ,  t  inn  arid  f  h  j  eshi  1 1  d  i  nq  ,  a  sincjle  photon  sometimes  results  in  more 
to. if-  i o e  pix.-L  r »-c o r *1  i nt)  a  one,  depending  upon  the  threshold  level  and  the  size  of  the 
t  I'.s  f  -  ♦  ligt.f  '-/itinq  from  the  image  i  n  t  eri  s  i  t  i  e  r  .  Tatile  1  shows  the  autocorrelations 
jos  Mo-  ir  •dividual  sipidred  transfer  tunctions  of  some  of  the  various  patterns  of  ones  re- 
i  j  *  i  f  -  ;  fr  nn  o  s  i  rii]  1  »•  pfmt'in.  Fart*  pattern  is,  in  effect,  t  f»e  impulse  response  of  the  de  - 
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Table  1.  Kvent  detection  data:  individual 
impulse  responses,  their  auto  correlations, 
and  their  power  spectra. 


DFTFCTION 

IMPU1.SK 

RKSPONSF 

AC  OF 
IMPUI.SF 
RESPONSE 

DHTFCTION 

TRANSFER 

FUNCTION' 

m 

□ 

1 

UH] 

E 

HH 

2  ♦  2co M2  *  u/N) 

0 

1 

2  +  2cos(2^  v/N) 

HU 

Bum 

3  ♦  4co& (2n  u/N) 

+  2cos(4  it  u/N) 

\ 

m 

3  +  2cos<2^r  u/N j 

llllj  m 

? 

ID 

♦  2cos(2Tr  v/N) 

0 

_L 

♦  2cos  ^2  *(u+v)  /N  j 

Figure  2.  Summed  power  spectrum  of  an  un¬ 
resolved  star  (linear  scale). 

The  middle  and  upper  curves  are 
the  same  as  the  lower  curve,  ex¬ 
cept  have  10X  and  100X  vertical 
scales,  respectively. 

te  Men  system;  and  in  any  one  image,  several  different  patterns  may  appear.  That  is,  this 
impulse  respur s *  may  vary  from  photon  to  photon  within  the  same  image.  Assuming  a  sparse 
population  of  photons  within  each  image,  it.  can  be  shown  that  the  net  squared  transfer 
function,  ,'je  ,  the  .  r.«,en>hie  of  pho t on-produr  *»d  patterns  within  an  image,  is  given  by  a 
weighted  sum  of  the  individual  squared  transfer  functions  of  the  individual  patterns.  We 
ref  “l  n  this  w>  iqht-'::  sum  as  the  deter  t  ion  transfer  function  squared  (DTF?). 

Om  can  er.npeny  ate  tor  the  noise  bias  term  hy  the  following  steps'?.  (1)  Over  the 
spa'  iul  fr'-qut-ni  ies  above  the  telescope  cut-off,  perform  a  two-dimensional  least-squares 
fit  t  i  i  qh  t  >■  ;  sum  of  individual  squared  transfer  functions  (some  of  which  are  shown  in 
Title  1)  t  fee  uiimh  ii  power  spectrum.  by  this,  the  DTF?  is  determined. ^  (2)  Compensate 
fi  i-f'-rt  .*  Inc  fiTf  2  py  dividing  the  summed  power  spectrum  by  the  OTF?  (for  all  spa- 

*  iai  frequencies).  Hy  this,  the  noise  bias  term  is  made  a  constant.  (?)  Subtract  from 
t  ’  :  !‘  ‘-rim,;.,.:  cummed  power  spectrum  the  constant  noise  bias  term.  This  DTF? 

. i r. i ;  niise  bias  compensation  are  demonstrated  in  Figures  3  and  4  for  the  binary  star  system 
’  41  In  this  case,  the  magnification  was  3UX  and  the  wavelength  was  750  nm  (10  nm 

spectral  bandwidth)  and  so  the  telescope  cut-off  is  at  a  spatial  frequency  of  74  pixels. 
This  data  set  resulted  from  power -spec t rum  averaging  of  1820  short  exposure  images  contain¬ 
ing  a  l  tal  of  about  2.4  x  105  photons. 


I 'i  tig-  aut  in  nr  r  el  a  t  i  on  domain,  the  noise  bias  term  results  in  a  spike  at  the  (0,  0) 
i.  ..ordinate,  and  the  OtF?  causes  the  spike  to  he  spread  over  a  few  pixels  about  (0,  0). 

.-in  pens*  t  ion  for  the  DTF?  causes  the  spike  to  collapse  to  a  de  1 1  a- f  unr  t  i  on  at  (0,  0). 

Then  the  subtraction  of  the  noise  bias  In  the  Fourier  domain  removes  the  del t a- f unr t ion  at 
(U,  U  /  iri  the  autocorrelation. 

More  generally,  the  functional  form  of  the  DTf ?  is  heavily  dependent  on  the  manner  In 
which  the  images  are  detected  and  should  he  modified  according  to  the  characteristics  of 
tin-  dete.  t  inn  hardware  used. 

2 

1  f  p  HI  k  ]  („•  M  T  F  C 


i.  fj(nppns  tion  for  thn  speckle  MTK?  would  ordinarily  he  accomplished  hy  dividino  the 
power  spectrum  hy  t ti e  summed  power  spectrum  of  a  reference  star!.  Both  power 
spef.fr  a  should  first  he  corrected  fur  the  DTF^  and  the  nnise  bias  term. 
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t  M.  *  t» 

PIXELS 

it  spectrum  of  the  binary  SAO 
t>  i .  ( A )  Middle  curve  raw  s urn - 

power  spectrum;  (B)  upper 
vt  summed  power  spectrum  di¬ 
ed  by  r.  he  DTK^ ,  (C)  lower 
vt*  imV  -  compensated  summed 
er  spectrum  with  noise  bias 
■  r acted . 


Figure  4.  (A) -(C)  Two  dimensional  view  of 

Figure  .3  (A)  -  (C)  ,  (D)  the  Fourier 
modulus,  i.e. ,  the  square  root 
of  (C) .  Note:  the  residual 
noise  beyond  the  telescope  cut¬ 
off  frequency  is  visible  in 
this  case  and  not  in  (C)  because 
the  square  root  operation  re¬ 
duces  the  dynamic  range  of  the 
da  in . 
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Figure  5.  Worden  subtract  method  on  an  unresolved  star  (low  spatial 

frequencies).  (A)  upper  curve:  summed  power  spectrum;  low¬ 
er  curve:  power  spectrum  of  the  sum  of  the  centroided  im¬ 
ages;  (B)  upper  curve:  summed  power  spectrum  (note  expanded 
vertical  scale),  lower  curve  summed  power  spectrum  minus 
the  power  spectrum  of  the  sum  of  the  centroided  images; 

(C)  same  as  (B) ,  except  a  smaller  percentage  of  the  power 
spectrum  of  the  sum  of  the  centroided  images  was  subtracted. 

"  :  •  '  i  i  h  would  (•>•  true  fur  objects  of  diameter  only  a  small  fraction  of  an  arc-sec), 

«•  re  |  1  <•  .!  power  -pe  -  t  rum  in  that  reel  ion  by  a  constant.  The  constant  i«  chosen 

t  !  •:  ■  •„-•••  i  ,  tent  with  the  value  <f  the  summed  power  spectrum  in  the  region  just  beyond  t'ne 
«•  j  v  -  *  -  f  r  i  uu-n -y  jjiwe.  flr  in  t  n  i  Case  of  the  Worden  subtract,  method,  this  method  due' 

•  rru't  ter  Ue  r.i  du  1  e  -  f  i  enuenr.y  vs.  h  i  qhe  r  -  f  r  equency  regions  of  the  speck  le  MTF  2; 

■  ..t  r  ,  as  rated  earlier,  t'e  >  perk  I  e  MTt  2  is  reasonably  well  behaved  for  those  spatial 
t  r  e,;  ,;pn  r  jes ,  .mo  correct  ing  hr  the  very  low  spatial  frequencies  corrects  for  the  greatest 
,.  j  :  •  f  t  'ii  i  r  r  oi  . 


:•  •  si  t  it  lift  in.)  the  summed  ;n  w-  r  spect  rum  to  curiect  fot  the  speckle  MTF 2  |s 

!•  •>  ur's  t  rated  in  f  inure  w  tor  the  ninary  '.AO  0416?  for  which  reference  star  data  was  not 
i.iiiatle.  I  '•  order  to  inert;  ase  the  arr  m  ar  y  of  the  assumption  tliat  the  Fourier  modulus 
( ...  r  it.-,  square,  the  power  spectrum)  is  constant  for  very  low  spatial  frequencies,  the  OTF 
a*  '  "i.i  ,c -h  i  as -1. -j  r  r  er  t  pi)  Fourier  modulus  was  divided  by  the  MTF  due  to  the  telescope  aper- 
t  j  r>  !  *n  i  r:h  was  approximated  by  the  MTF  due  to  3  circular  aperture  of  diameter  2.3  meters). 
I-'  1  llipti'  al  shape  of  t  he  Ic.urier  modulus  data  is  due  tn  the  difference  in  scale  factors 

1-  the  two  dimensions  as  nuten  earlier.  within  the  low  frequency  region,  wherever  the 
i  ..tier  modi;  I  us  ext  r.t’ded  a  threshold  value,  it  was  clipped  tn  that  threshold  value.  The 
re.ult  was  multiplied  by  the  MTF  due  to  the  telescope  t.o  arrive  at  our  final  estimate  of 
t'»'  f  '11 1  r  i  ■  modulus  of  aAO  ’Min'.  including  the  telescope  MTF.  In  the  process  of  multiply¬ 
ing  hack  in  the  telescope  mtf,  the  residual  noise  beyond  the  telescope  cut-off  frequency 
was  set  to  zero. 

image  reconstruction  results 

The  Fourier  modulus  estimate  shown  in  Figure  6(d)  was  truncated  to  a  128  x  128  array, 
in  or  ter  to  save  computation  time  in  the  iterative  reconstruction.  This  caused  a  slight 
truncation  of  the  highest  spatial  frequencies  along  the  horizontal  dimension  of  Fig- 
j  1  e  6(d).  s  A  ( 1  *4163  was  reconstructed  ijs  Inc  The  iterative  method,  and  the  imagi  s  resulting 
from  two  different  selections  of  toe  initial  input  tn  the  algorithm  are  shown  in  Fig¬ 
ures  7(a)  and  (h),  respectively.  The  rms  error  F0  was  reduced  to  about  O.OS.  For  the 
purpose  of  display,  a  (sin  x)/x  interpolation  was  performed  on  the  images  of  Figure  7  in 
urlei  to  increase  the  sampling  rate  across  the  image.  In  order  tn  get  an  indication  of 
the  sensitivity  of  the  method  to  the  clipping  threshold  level  described  in  the  previous 
"i  tion,  toe  c  lipping  was  done  over  again  using  a  3 1%  greater  threshold  value  (which  is 
obviously  greater  than  the  optimum  threshold).  Two  images  reconstructed  from  the  result inq 
Fourier  modulus  estimate  are  shown  in  Figure  7(c)  and  (d).  Half  the  time,  the  iterative 
r ei  (instruction  algorithm  produces  an  image  rotated  by  180°  due  to  the  inherent  2-fnld  am¬ 
biguity  uf  the  Fourier  modulus  data. 
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Fir.uro  h.  C  l  i  pp  im;  to  compensate  lor  the 
speckle  MIT  (for  "see i no")  for 
the  binary  SAO  (A)  Four¬ 

ier  modulus,  same  as  Figure  ad)), 
{ H)  Fourier  modulus  compensated 
tor  t  e lescop e  MT F  ( a 1 1  etnp t. e d 
division  by  ::ero  is  evident  tor 
spatial  frequencies  above  the 
telescope  cut -oft  frequency). 

(C)  clipping  of  the  low  spatial 
frequencies;  (D)  Fourier  modulus 
estimate  obtained  by  putting  back 
in  the  telescope  MT F . 


Fiy.urt-  /.  Reconstructed  images  of  SAO 
R 'd  M  (  see  t  ext )  . 
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••  •  s.io-  :  i  in  t  n  raw  sun  rif.-it  po#n-r  ■  p'-stium:  mdi-  !«•••»  truv.f- 1  tor  f  ;  i  '  s  i  •  -mrer  ■  \  t  i  u« , , 

•  i  '•»  tie  .  ,  t  •  t  :  i  i  i 1  m  ,  - » r ;  I  • .  p  ►  •  r  *  J » •  f  *  7  f-  rrri.p*  ns  a  *  j  ef  .  f  <  r  ■  '•fr ;  •  r  •-  *  i  •  .  t  f  •  •  •  v  1 1 

rf't'.renr*-  star  Oat  j  is  not.  available,  ari  improvement  r.v*»r  the  Aorc*- r.  <- u- 1  r  •,  r  t  nmt.hort 
t '  *-  •  i"ple  fn-thro:  f  »- 1 1  pp  i  mj  *  *'#••  t  «•  o  r  ie  r  ■  r?  i : ! « ;  ■  •;  1  >  *  i ;  ’*•  v'  :  v  i  •  w  pit  id.  ‘rrT.,rr- 
i  'for  of!  jet  t s  much  smaller  than  one  arc-sen  in  iti.ifnrt«r ) .  A  ur  st  nn't  i  on  nf  the 
■  i  •  or.  i*.'*  u'driu  t  \  is  r.»  th-rj  r^-sul  t^-ri  in  nr  im.-nn-  r.,v  i  r-u  t  t  i'nr,  •  '-eat  d  i  ■ .  if  f  7 

a  r  - '  t  it  an  unde  of  uS./0  (47.3°)  arnt  a  magnitude  o  i  f  t  •  r  r-:  u  »•>  .d  Id-7  (brightness  ratio 
f  Tre  r  er  ore,  r  ruf  t,  i  on  of  a  tjjnary  i  •'  trivial  .me  :*u.sir-»  tn,  .t  t^e 

iterative  ’n^thorj;  however,  the  iterative  method  is  required  t^r  <  •  mp  1  \  rot  >-c  •.  ty'ec  t  s ,  and 
‘  '  •  :-d  !  ;r  .  'do-p.s  dr- <  r  i  h»*ii  t.ere  t>u  t.  hi--.  • .  i  r»-  r  d  •  >  irq  e  :*  »  4  d-  u*  'r>-  om- 


i  r  r  I.ifird  nnr.  OS  .  Hosed  on  the  ojrr  p*/;  .it  imad*-  I>o  nn  d  ru-'  f  i  on  i  x  |  •  r  i  m*‘ u  t  i  m]  .  jmm.>  d 

;  f  ;  t  r  ->  - 1  r  ■  ■ . mp  ]  i  r  ,,  * ,  <\  t  wi  -cim'-iv  inn.-d  oh  .  ?  >  i  ry  i .  *  •  :  -  i  m >  i  !  \  t  ■  .  ■  f  i  ■  ■  1  .  v  t  ‘ 

i  T  ;  j  r  -  *  •  r  i  t  srbul  er.re  and  |  t-r  t  ,  ;f .  no  j  .•  ^  ,  i  ‘  i  •  * ;  •  i  ?  ■  ;  *  •  •  *  i  *  will  be  ;  e  - 

t  r  »■:  ons  1 1  us  t  *  i  n»  -r»  so  1  of  \  on  imacjes  uf  roinpl  in.it  >  -4  t  t  r  ■  m  i  <■  a  ]  » ■  t  ■  . 
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Abstract 


're  uniqueness  of  solutions  to  the  phase  retrieval  problem  is 
explored  using  the  theory  of  analytic  functions.  It  is  shown  that 
:r  a  function  or  its  autocorrel ation  satisfy  certain  disconnection 
conditions,  then  the  solution  is  unambiguous  unless  the  separated 
parts  of  the  function  are  related  to  one  another  in  a  special  way. 
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<  .  ;  ntr  Auction 

"o~  pr .  u  i  e:n  of  phase  retrieval  arises  in  many  tu-lds;  s  o»- 1  a  1 
is  'r'lvn  v ,  racio  astronomy,  radar,  antenna  theory,  holography,  inter 
* e r . ant t r y ,  crystal  loqrapny,  electron  microscopy,  conr  ier  transt-rm 
spe ..  t  r  js  i.  o  py ,  and  snape  recognition  'Fourier  descri ptors ) .  o  yen  a 
■  1  . 1 1 • ' n  to  the  phase  retrieval  protlem,  one  would  i  ue  to  mow  it 

•  t  1  s  <t  •  o  i,e .  *ne  general  one -dine  os  ional  phase  retrieval  i,niquenes 
pr  ';,i  if  js  jtateo  as  follows:  ’ f  f  is  a  rompl- <-va  1  ,eci  function 
■VT  ’neo  on  tne  real  Iv  e  and  F  is  't.s  Fourier  transform,  under  what 
o.n.rt'^ns  aid  to  wnat  extent  >s  f  ietermined  py  p  ?  Since  F  is 
t  n  -  F  u r  '  er  tries  form  of  the  autocorrelation  fur-cfor  at  r,  -his 

or  ;z  '■  em  ;  .r:  :  e  jqui  ;a!ently  stated  in  tne  form:  'aider  what  conn  - 
and  to  wn.it  extent  ’S  *'  determined  py  its  autocorrelation 

•  .not’  on  ? 

i  r.  th’s  PrfD'xr,  restrict  Pur  attention  tc  the  i.ase  in  which  f 
•s  assumed  to  nave  compact  support  (i.e.,  f (xl  *  0  tor  <  outs  ide  or 
some  *imte  interval).  Metncds  of  reducina  the  nuinber  of  tjnct’mi 
r  na.-  ing  tne  same  f  by  imposing  disconnectedness  conditions  or-  Sr> 
s,.  r.ocrt  at  ♦  ;r  toe  support  of  its  autocorrelation  *  .met  ’  on  are 
::"jren.  S^'oe  .■  r  tne  theorems  presented  in  this  paper  are  refinement 
-.i  :  ;r  :er.pr  ai  i  z3tions  of  xnown  results  and  sore  are  n*w.  ->!!  proof 
ar--  "•  tne  appendi « . 

.  '--rmmO  M.  qy  ana  0<-o]  imirar ’es 

’r.e  support,  Z-e%  of  a  complex-valued  measurable  (  .r.ct’cn,  or. 
,'*iai  ''ne,  H,  1  s  tne  smallest  closed  subset  ,r  ».)e  real  ’ire 
•  t'.’  >  *  wn  ■ : >:  *  is  zero  almost  everywhere  (a.e.l.  tince  tv  •: n an o 
•nq  -t.s  x  -.lues  on  a  set  of  measure  zero  a  function  can  be  mace  equal 

•  i  z-r o  everywhere  oft  its  support,  in  the  remainder  of  tn is  paper 

iil  *  uri  •'  *  i  or.s  re  assumed  to  be  zerr>  everywhere  off  tneir  sup¬ 
ports.  'ne  'titer',  i.  of  Support,  I  .  f  =  '  af ,  bf  1  ,  or  ♦  is  the  jrra  1  i  - 

••St  '  !  .  s  ed  :  1 1  *>■  >'  .  i  ;  ,r,  f  *  i  "  i  n  q  j,  .  ”ne  er'  te  r  ,  r  »  t  '  ;  ]  j  <-ri 
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:  ,,r  >,et  of  H  ">  compact  1  f  1 1  1,3  contained  in  some 

%..•■.  j|_  C  's  tne  crimp  lex  plane.  If  wr  if,  then  w  is  its 
1  i  is  the  space  of  all  complex-valued 

,nr  an  i  e  •  ,n  t  ions  on  H  with  compact  support. 

-  \{  i,  the  involution  of  f  ann  q  is  given  by 

*  *  1  <  ■  =  f  '  /  j  q  (  x  -  y  I  dy . 

*  t «o  t.jnrt  10ns  m  ~  1  fi  )  is  also  in  '_0(  5  ).  We 

.  -  ♦  . .  .  rne  au tocorre  la t ton  of  f  is  given  by 

iu  to  ( f '  =  f  *  f . 

.  •  ,•  r  .x*->n  xf  *  >s  li  ven  b  y 

-  w  =  f  <  e  ' W  *  px  <£. 

.^tt-vs  „ • ' 1  De  .sec  to  denote  tne  Laplace  transforms  of  the 
■>in  r./  tne  '.or responding  lower  case  letters.  If  F  is 
.■  ,r  if  :»  > ^  le'OH]  oy  F*iw)  =  F(w).  If  F  is  the 
•  t  r-.  *  ♦  t  ->irn  c*  is  tne  uaplace  transform  of  f.  Also, 

U'.;J  ♦  r.  .  *  -r<~-  t  0,  '■hen  the  Laplace  transform  of 
.  .  ,  *ne  _  -  p  ljce  transform  of  auto(f)  =  f  *  f 


Mri'sn  jn  (f  ,  tne  function  -c  's  defined 


-  \£ 


.•..real  ter,.'  0*  F  tnen  nc 1  w  i  'S  ’t  or  ter. 

>t>r  ■)  ijt  £  tnen  npiw  i  'S  one-nai*  ’  t  $  Orcer , 


K  •  , 


1 


1 


'V 


0(F 


j  w  .  mu";  >  u ,  npp 
"piwl  =  nc(w)  -  rip(w). 


w:  Im(w)  >  0,  nPp*(w)  >  0  and  w  *  0  (■, 

F< 


an 


Jp(W)  > 


; f  A  and  B  are  sets, 
f 

a  '  3  =  jx:  <tA  ana  x<B 
:f  a,  B  =  H,  then 


-  s 


=  |  a  -  b:  ae.A  and  beB 


tre  complement  of  A  in  K  . 


!t  a  and  c  are  real  numbers,  the  functions  f(x),  f(x  +  a)e1c  and 
-'.x  *  aie1c  all  nave  the  same  Fourier  .modulus.  If  these  are  the  only 
'.nations  with  that  Fourier  modulus,  then  f(x)  is  said  to  be  un  ique 
jr  i  ’ts  Fourier  modulus  is  said  to  be  unambiguous .  Otherwise,  f  is 
"i.)n -unique  and  its  Fourier  modulus  is  ambiguous .  If 

g ( x )  =  f ( x  +  a)e’c  or  g ( x )  =  f(x  +  a)eic 
‘  ■  en  f  ;nd  q  are  equivalent  or  in  symbols,  f  ~  g. 


Examples  at'  *ion-ijn iqueness 

.et  be  a  positive  integer  and  A  =  ^1,  ?,  .  .  . ,  N^-.  Let  B  be 
a  ..  reset  of  A  and  C  =  A  \  B. 


’ueorem  Let  b  n  =  1,  .  .  .,  N,  and  a  be  given  complex 

,rper s  and  let  c  and  d  ,  n  =  0,  1,  .  .  .,  N,  be  defined  bv 


■n< 

n  =  1 


N 

V"  n 

/  .  cnx 

n=0 
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in  j 


uet 


Ui-K*  * 


ncB 


1) 


n 

nEC 


( x  -  b  ) 
n ' 


an  arD-'trary  function  in  L^(  TJ 

o 


N 

VV  n 

_ _ .  d  x 

n=0  n 


,  8  a  real  number , 


qi 


V  .  (l 

_  V  'x  -  nB  ) 

n=0 


in 


n=L 


►-  f 


.  G  i  u  1  =  u .  .1  V-  ije  J\  . 

I 

'•  t-:  a  sellar  r  actor  i  za  t  ion  rero.Maue  ’s  used  Dy  8rUC<  anc 
|  - x  imp !e  1 :  bet 

1  for  x  <  1/2 
0  otherwise 

A-r  ■’  ;v- 

1  *  3)  (x  +  2)  =  6  *  5x  +  x? 


a  1  =  or'  x  +  5 f  ■:  x  -  &\  *  f  ( x  _  s) 


Sod’n 
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n  (  x  )  =  2 1  ( x  )  +  7  f  (  x  -  4 )  +•  3  f  ( x  -  3  ) 

T  hen ,  by  "heorem  1,  .  G  (  li  )  ]  =  ,H(u)j  ¥  ue  (see  Figure  1). 

(n  tnis  example,  the  two  functions  are  quite  similar  in  appear- 
iixi.  In  particular,  they  nave  the  same  support.  In  a  particular 
application,  it  would  probably  oe  difficult  to  rule  out  either  one 
a-,  "ne  r 1  gnt  solution  on  the  basis  of  a  priori  information. 

: '.ample  2:  _et  f  oe  as  in  example  1.  now  we  use 

2  <  l '  ( x c  -  i  x  +  2)  -  8  +  2 x  ♦  3 x'~  2 x 

2 !  2  x  *  1  ; !  x -  --j-x  +  2)  =  4  +  7»  *  4x'J. 

m  <  =  ofix)  +  2 f ( x  -  4)  +  3f / x  -  3)  +  2f'x  -  12) 

:f  >1 

n i x  )  =  4  f ( x  )  +  7  f  f  x  -  41  +  4  f ( x  -  12) 

--  :  ’  o. ..  x- ~  .  Pca’ii,  oy  theorem  1,  6(u)  =  M(u)'  ¥  Ue  ft  . 

e>ampie  snows  tnat,  jr  1  ike  the  situation  in  Example  1,  the 
•  *  -  ->ns  ;  an  nave  different  supports. 

■■  m  * /,  1  xinq  theorem  gives  another  method  for  generating  exam- 
_  -rs  f  '■  )'i- uniqueness. 

’’•eur-m  2:  .et  *,  qeL‘  (  ft  )  and 

h,  =  f  *  q 
n.  =  f  *  g. 

-  <  U  =  hi...  ■  UJS-  ¥  Jr  1\  • 
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'V 


examples  1  and  2  could  have  been  obtained  by  using  Theorem  2.  In 
gceril,  ’’heorem  c  is  useful  in  generating  examples  of  non-uniqueness 
’n  .vnicn  oot.n  functions  are  positive  because  if  f  and  a  are  positive, 
teen  ,o  art  n.  and  n;.  (bote:  Theorem  2  also  holds  for  functions 
or  r wo  or  more  variables.) 


urr'eS  or  uniqueness 

;  jra  T  are  presented  in  this  section  out  of  logical 
’•  :.r-  :  r-r j  are  needed  in  t.ne  discussion  of  examples  of  unique 

:■  ■  .  pruofs  l  see  appendix)  use  Corollary  3  to  "heorem  3 

:  to  "neorem  3,  a  function  t  is  unique  if  r  r,as  no 


£  ■.ample  3:  _et  f  oe  as  in  Examples  1  and  2.  Then 

It  s i n  w  n 

2  -  for  w  p  0 

w 

2  for  w  =  0. 

:  ’  is  no  non-real  zeroes  and  therefore  f  is  unique. 

"-e  -'c-i  lowing  theorem  gives  a  method  for  generating  more  examples 

,n  ’ q  jeress  . 


"  ec  r  em 
•  s  a  i "  n  : 


.et  fcL  (  33  )  be  unique  and  let  bn  be  complex  num- 

1,  n  =  1,  .  .  .,  b.  Let  cn ,  n  =  0,  1,  .  .  .,  b,  be 


n-  -  bn ■  = 


n  =  l 


n=0 


an  j  i-t 
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q(x  ;  =  Y"  C  f ( x  -  en ) 
n=0 

where  t  >s  d  real  number.  Then  q  is  unique. 

Cord  i  ary:  If  f  is  unique  and 

N 

g ( x )  =  -  en) 

n=0 


then  c  is  unique. 


Lx  amp  1c  4:  Let  f  be  as  in  the  previous  examples  and  let 

IK 

g(x)  =  £f(x  _  4n)- 
n=0 

Tner, ,  d>  tne  corollary  to  Theorem  3,  g  is  unique  (see  Figure  3). 

Another  method  for  generating  examples  of  uniqueness  is  given  Dy 
tne  following  theorem. 

2 

Theorem  4;  If  f,  ge L ^ ( 3S  ),  f  is  unique  and  6  has  no  non- 
real  zeroes,  then  f  *  g  is  unique. 

Example  4:  Let  f  be  as  in  the  previous  examples  and  let  o  =  f  r  f. 
~ner.,  Dy  Theorem  4,  a  is  unique  (see  Figure  4). 

b .  Factor  i  zatior.  of  F 

Let  fcL^(  K  ).  Then  auto(f)  =  f  *  f  e  L4(  )  and  it  follows 

from  the  Paley-weiner  Theorem  ([2],  p.  103)  that  both  F  and  FF*  are 
entire  functions  of  exponential  type.  Since  FF*  is  entire,  D(F)  (see 
Section  2)  is  countable.  Let  its  elements  be  numbered 
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0(F)  = 


w  :n  =  1,  ?, 
n  ’  ’ 


Such  that 


U' 


is  a  non-aecreas ing  sequence.  The  following 


factorization  is  essentially  due  to  Titchmarsh  ([3],  p.  285,  Theorem 
Vi).  (Note:  In  [3],  Titchmarsh  defines  the  Laplace  transform  as 


F  (w)  = 


r  ,  V  WX  „ 

f ( x )  e  dx 


as  opposed  to  our 

F(w)  =  f  f(x)  e  x  dx. 

-/—DC 


Tnerefore,  in  using  his  results,  appropriate  adjustments  must  be 
made .  ) 


Theorem  5:  we  have 


(w)  =  Ij-  F(m)(0)w,M  e 


■  "CfW  n 


n  =  l 


nr(wr 


1  -  — 


nF(V 


where  rr  =  2  n  p  ( 0 )  ana  the  infinite  product  is  conditionally 
convergent . 


Tnis  particular  form  of  factor ization  of  F  is  chosen  to  facili¬ 
tate  "zero  flipping"  arguments  which  will  be  used  later. 

in  the  sequel,  we  will  also  need  the  following  theorem. 

Theorem  6:  If  6  is  an  integer-valued  function  on  (T  3  0  <  B(w)  <  r 
¥  wc  €  N  Band  0  <_  6  ( u )  <_  2np(u)  ¥  ue  B  ,  then  tne  infinite  product 
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6  'wn 


n  — ^ 

i  i  -  — 
n=l  w 

1—  n 


is  absolutely  convergent  ¥  we  (C  . 

Tne  next  theorem  has  to  do  witn  "flipping  zeroes."  More  v. i  1 1  oe 
said  about  this  in  the  next  section. 


T heorem  7 :  Let  feL^  33  )  and 


f  N 

i  w 

n  , 
i  Jn= 


be  the  distinct  non-real 


'  'IIS  I 

zeroes  of  F.  Let  8  be  an  integer- valued  function  on  £  3  0  <  e(w)  <  r,  (w 
for  wt  (f  \  -K  and  0  <  6  ( u )  <  2r,p(u)  for  ue  33  .  Let 


G  ( w )  =  F (w) 


1  -  *L' 


1  -  - 


B(wn) 


anc  let  g  be  the  inverse  transform  of  G.  Then  ;G(u)’  =  ,F(u)  for 


uc  33 ,  geL  (  1a  )  an  a 


Kg)  =  l(f). 


6 .  Functions  with  Disconnected  Support 

In  tnis  section,  it  is  shown  that  the  probability  of  uniqueness 
is  much  higher  for  functions  whose  supports  satisfy  certain  discon¬ 
nectedness  conditions. 

Let  I  n  =  1 ,  .  .  . ,  N,  be  a  sequence  of  disjoint  closed  inter¬ 
vals  satisfying 

6.D  dn  -  V  n  (ij  -  ik)  -t 

for  j  4  f  and  <n,  m>  f  <j,  k>  (where  <  ,  >  denotes  ordered  pair). 

Let 
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A  =  U  I. 

n 


-et  f ,  QeLp(  1\  )  satisfy 


anc  let 


S(f  )  c  A  and  S(g)  C 


V*> 


f(x)  for  xcl 


0  otherwise 


f  or  n  =  1 , 


(x)  for  xt I 


9n(x) 


^0  otherwise 


N.  Finally,  let 


B  =  O  2 (Fr 


assume  here  that  f  ^  0,  n  =  1 ,  .  .  . ,  n. 


Theorem  8:  Let  f,  g,  f  and  g  be  as  described  above.  |F| 


lG(uii  V  uc  B  iff  3  a  real  number  e  and  an  integer-valued  funct1 

a  defined  on  (C  with  0  <  a(z)  <  min  [n  (z)]  f  zc  £  \  ft  ana 

l<n<N  rn 

C  <  a ( u )  <  2  min  [rtp(u)]  for  uc  33  such  that  if 
l<n<N  n 


u>  = 


o(w)  =  e 


1  -  ^ 
2EBL  z. 
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then  tor ‘  f.'  /  ?, 


i(Cf-C  )w 

1  )  Gp(wl  =  e  9  <*'J(w)fn(w) 

ana  cf  -  C  *  cf  •  c  for  n  .  1 .  N, 

n  yn  y 

ana  for  n  =  2,  either  1)  holas  or 

*  i(cf+c  )w 

..  !  o(  U  :  =  e  -  C>(w)F  (w) 

an cJ  c ,  *  c  =  cf  +  c  for  n  =  1,  2. 

f  q  f  q  ’ 

if  the  'unction  f  is  a  more  or  less  arbitrary  function  gotten 
from  trie  real  world  (but  satisfying  the  hypothesis  of  Theorem  8), 
tner  it  should  almost  always  be  the  case  that  E  =  t.  For  this  case, 
we  obtair  tnt  following  result: 

,C"\r  lai'>  i:  it  E  =  t,  then  'F(u)  =  i G ( u ) .  V-  ut  K  iff  f  '  o. 

In  [ij,  Bates  states  a  similar  result  but  with  too  weak  a  hypoth¬ 
esis  ar, a  too  strync  a  conclusion.  (See  both  the  discussion  lmme- 
ciate'v  following  ano  the  discussion  following  Corollary  3  to 
Theorem.  S.)  it  should  be  notec  that  Corollary  1  does  not  hold  if  it 
is  merely  assunec  tnat  the  I  are  disjoint.  To  see  this,  let  f,  c 
ana  n  oe  deemed  as  in  Example  1.  Let 

=  [-0.5,  0.5],  1 ^  =  L^.5,  5.5]  and  =  [7.5,  8.5]. 

Let  q  ana  h  play  trie  roles  of  f  and  q  in  Corollary  1.  Then  Z(G  )  =  t 
f or  ri  =  1 ,  ?,  3  anc  hence  B  =  i.  Also,  ;G(u)  =  H(u)  ¥  ue  B  but 
u  ana  h  are  not  equivalent. 

By  setting  N  =  1  in  Theorem  8,  we  obtain  the  basic  "zpro- 
f lipping"  result  of  Hofstetter  and  Walther. 
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:  or  o  1  1  a  ry  2'  (Hofstetter  [5j  -  waltner  [  P  j ) :  F  ( u  )  =  G  ( u  1 

V-  ut  ft  i  ff  3  a  real  number  e  and  an  integer-valued  function  a  de¬ 
fined  on  Z(F)  with  0  <  a(z)  <  r,  ( z )  ¥  ZcZ(F)  such  tnat 


i(c,-c  )w  TT 

>1&  e  f  9  Ffw)  I  1 

ZcZ(F) 


Flipping  all  the  non-real  zeroes  of  F  yields  tne  following 
relations : 


Since  F*  is  tne  Laplace  transform  of  f  and  f  is  equivalent  to  f, 
we  obtain: 

Corollary  3:  Tne  followinq  are  equivalent: 

1.  F(u)  is  unambiguous. 

2.  F(w)  has  no  non-real  zeroes  or  one  non-real  zero  of  order  1. 

3.  F(w)F*(w)  has  no  non-real  zeroes  or  two  non-real  zeroes  of 
order  1 . 
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oi"  i' 1 1  ary  1  suftprs  from  tne  ar  awDack  that  both  f  ar.c  a  must  lie 

N 

assumed  tc  have  their  supports  contained  in  A  =  ;  For  tnis 

n  =  ! 

reason,  't  cannot  be  concluded  that  f  is  unique.  Tne  necessity  cf 
t1  ’S  cond’tmn  car.  be  seen  from  the  following.  Let  f,  g  and  h  be  as 
de*  Hit"  in  :>ample  2.  Let  ri  play  the  role  of  trie  function  f  in 

1 "  e i 1  ary  ^  ■  -  1  So  let 

I  I  =  ,  -0.6,  4.1'  ana  l-  =  ,  1!.6,  12.6] . 

!  ,  arid  satisfy  the  separate'  confiitior  f.l.  Let  h  be  the 
ic  t  nir  of  n  to  I  ,  n  =  1,  ’her.  ,  _  (,  and  therefore 
:  -  f .  -is;  H(u)  =  i  G  ( u )  V-  uf  E  t-.it  i  aid  c  are  not  equivalent. 
'•  •  e ,  i  's  net  unique.  A  similar  result  occurs  ’n  tne  t  =  3 

i 1  •  *  *  a*'  .  •-  * 

-  i,  =  [3.5,  4.5]  and  I,  =  [11.5,  12.5], 


A*.'  t*  —  qraw;  a  •  that  tne  functions  f  and  g  of  Corollary  1 

're**  ‘yr*-  t*‘e  M.ow  leciae  that  their  supports  satisfy 
•  ;  •  •  •  .  ru..  *  on*  ii  tne  f  >rm  e*  a  priori  knowledge.  The  next 

:*  -  •"  :rv  >•:  *i  e  ■  >  or awt’dCKs  for  positive  f  by  imposing  a  separa- 
•  i'  t  ’  •  .*  '*  *-  support  of  the  autocorrelation  function  of  f 

,  ...  i , , t-  s i. pen  t  of  f.  Since  tne  autocorrelation  function  is 

.  a'  "  f *.*-  amt  for  any  two  functions  with  tne  same  Fourier 
■  tii*'.'-  v dif.pa c k s  disappear. 

:  -•  c  i ,*.rc  twr  lemmas.  Although  the  first  of  these  lemmas  is 

*  1»*  ,  *ur  positive  functions,  we  state  and  prove  it  (see 

apper1  c i »  f  or  art  itrary  conip  1  ex- val ued  functions. 


Lemma  2:  Let  feL^  (  E  ),  1(6)  =  [af,  bf]  and  I(auto(f))  = 

_  -  I- ,  cij .  Then 


d  =  bf-  af. 


66 


The  next  lemma  is  really  the  heart  of  the  matter. 

Lemma  3:  Let  tcL^(  B  ),  f  >0  and  1 ( au to ( f ) )  =  [-d,  dj.  Jf 
y  >  0  and 

r .  l  !  S  (  au  to  ('*’))  C  L~d,  -  ^  d]  U  (-y  ,  y)  U  l^t  d,  d  j 

tnei 

S(f)  C  [af,  a+  +  y)  m  (Df  -  y,  bf]. 

Note :  Lemma  3  is  not  true  if  the  condition  that  f  >  0  is 
dropped,  'ne  function  described  in  Figure  5  is  a  counterexample  for 

>  -  L  . 

'low  suppose  f  satisfies  the  hypothesis  of  Lemma  3  with  y  =  d/3. 


f, (x  )  = 


a'  i  c 


f  (x  )  for  xc[af ,  a ^  +  j  d) 
. 0  otherwise 

f ( x )  for  x e ( bf-  i  d ,  b 


f?(x)  = 


.0  otnerwise 
Using  tf!  i s  notation,  we  have  the  following  theorem. 


It 


arid 


Theorem  9:  Suppose  fcL^l  I\  ),  f  >  0  and  I ( auto ( f ) )  =  [-d,  d], 


S  ( au  to  ( f ) )  ~  [-Q,  -  d]  U  (-  j  d,  j  d)  U  d,  d] 


Z f F  -j )  H  Z(F2)  =  t. 
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t '  i  er  f  is  unique  a  mono  nonnegat've  functions; 
g  ^  G  anq 


G  ;■  u  )  =  F(u;  V-  u  €  I\ 


trier  q 


o 


;it  course,  us  me  ti'is  tneoreir.  requires  the  a  priori  Knowleaqe 
t > ■  t  fie  function  being  searened  for  is  positive.  However,  ther* 
r t  mar  y  applications  ii.  whirn  this  is  given. 


We  now  turn  our  attention  to  the  case  in  wich  G  is  gotten  from 

(  1  I- 


F  o v  f  iippinq  only  a  finite  number  of  non-real  zeroes.  : et  U  r 

1  njn  =  l 

oe  ai  arbitrary  *mite  sequence  of  distinct  non-real  zeroes  of  F  and 


Giv;  =  F  (w )  J~J 


n=  1 


w 

w 


Jr 


1  _  rL 

<—  hn  — 1 


whe'-e  trie  tr,  are  integers  satisfying  0  <  en  <  nF  f  w  )  ano  F  is  tne 

..aplace  transform  of  a  complex-valued  function  frL^f  I\  ).  et  q  be 

o  - ,  • 

tne  inverse  transform  of  G.  Then,  by  Theorem  7,  atLM  Ia  )  and 

;  i  f :  =  ;  <  a , . 

.et  '  ' *  •  ■  oenote  the  complement  of  S(f)  with  respect  tc  1a  .  Let 


S(f  )’  =  (-  «,  af )  L 


K 

U  p 

n 

m=l 


U  fcf,-) 


be  tne  aecompo^ ' t.  ion  of  the  open  set  S  ( f  )  ‘  into  disjoint  open  inter¬ 
vals  where  V  is  either  a  finite  inteqer  or  K  =  »  .  Similarly,  let 

"i  n 

Sfg  af)  u| 


J  0  I  U  (b.,  ■  ■ 


Us  =  l  J 
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k  >*  ’  i'r  ■  q  '  ■  >  ><  t-'t'.y  c  »  •> 

••  •  •  -  '■  ?’  .it ,  o .  ",  i  -  3  ar. ■  -  r  , 

c;  '  c 


t  fT  _  1  V 

P1  ’  V  ’  -  . . 


■Si.  unit  the  Clip  1 1  OPS  aetlnec  ’r  p  t-  pi- 


*  ,ir  s.  Pii  rr  3r  r 


y  t ,  ‘  pf  , 


Q  ,  ( 


j  f  ■  >  1  f  or  c,  <  y  < 

f  -  —  rr. 


otherwi se 


i"f !  >  )  for  t  <  >  <  Ljf 
i  rr,  -  —  f 


[  C  otherwise 


q ( > 1  for  a,  <  x  <  s 
t  —  —  rr 


otherwise 


c , 


fc,  (x  )  for  t(i  <_  x  <_  tf 


1 

lt>  otherwise 


trfcr  toe  wf  ,  n  =  1,  .  .  r,,  are  zeroes  of  both  f  ^  and 

r"‘l>-r  ■  _1  ar,« 


G  ( v. )  =  c  ^ 
r  r 


■n 

n=l  ! 


-  -  ^ 


i 1  - 


n— ‘ 


r  =  -  I  . 
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,  I  *  wa».  proved  t  v  ifee>.3»kt5  >•  „ 


toe  "  jS-  1'  vP  ,r' 


31:  C 


on  ones 


*\  -■ . .:  w''p'  a, 


. *..  *  X  ai’T  I"  V:':  " 

•  ’  te  r.  u'’". 1 lJ'  '  ’  art  f  1  '  PP'-'J  , 


v,t  ,  Wf  v  .Mis  idt  fu-'.-t  •.;»<-  l\  r w  •;>  *  r  '?«■ 

>0'  f,n-  lie*  ir-  )  t  1  C*f'  '■  ’S  3  *"'t1‘  S‘H  .  .  ■  P<srt  'C.i.d*  ,  fr  1  -  •- 1,11  ’• 1  ■ 

satisfied  !'  ‘  •  .IS  'H  ■  v  a  r.r-.’te  '  •'•:*■•  o<  Z‘-''c>es. 

Mu>  til  at  the  supports.  >'  ‘.urn  ♦  un :  t 1  ons  an-  tnr  *  nr' ' ' r’" 
s.i  1 1  srv  serf-  spent c  ’  t  • 


31  -  C 


:  o  r  S  l  :ri  p  i  1  C  1 


tv  ot  c-  scuss  ’or- ,  it  *  '  ’  t>e  ass  .liner  v  at  ■ 


'm  -i>er  intervals  '°r  =  (sm,  t^'',  *  =  '  , 
a:.;  o  v  f  . 


i> ■  t  1  r,p: 


er 


Assume  Cf  =  0.  :f  v.'c  ’•  a  finite  set  an, 


i  f  I  x  )  for  >  <r  0 
f  ,  •'  x  )  =  •: 

!  0  otne-wi se 
1  f • . i  for  x  >  0 

=  i 

;  o  otherwise 

1:  a  zero  of  tiotn  F  1  ana  F]  of  orders  > 


-  - ;  z 1  is  t 


\  ’  ,r  *  r  Of*  in  It  1  or.  of  Cp  ) 


ana 


i'p. 


(  r 


P  ,  U.  =  — j 


i,w  n 

F""°)  z  t  W  ( F  i 


rl  -  “| 

~5j 


=  r‘  r  (  0  . 
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'neorem  K  :  Assume  =  fl.  If  k(F  )  is  finite  arc  for  some 


>  0  ana  P  Pi  [-P  i  /  t  then; 
m,  m. 


l^f  ( x  )  for  x  £  -t 


m,. 

z 


f  i  (  x  )  = 


V’>  =  1 


0  otherwise 


f (x )  for  -s  <  x  <  S 

rrip  -  —  rru 


_0  otherwise 


f ( x  )  for  x  >  t 


f  ,  (  x  )  =  •< 


1^0  otherwise 

■en  each  ze*(F  )  is  a  zer -  of  F_ 1 ,  Fo  and  F 1  of  orders  >  Cp ( z ) , 


-  •.  tv.; 


F  '  r  1  I  0  i 


^  F,w  n 

ZeW(F) 


'-fl 

_ z  ; 

1  -  * 


oF(z) 


an 


F*(wi  =  F  ( 

0  F 1 rT7o)  0 


w)  n 


W~1 


0  F  '  z 


!  i  _  * 
zcW(F )  ]  z 


where  r  =  ?rip  (  0  )  . 
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Mi  mm  a  r  v 


Tr<t  pr  v'lilt't:  j*  jn  iqueness  in  phase  retrieval  was  ~yp ior*c . 
n a,:;.  <>■  ueneratinq  unique  and  non-unique  functions  were  p 
was  sh'.iwn  that  if  a  function  is  comprised  o*  parts  w’tr  in  real  or ‘ 
suppo*'t  Sut  sfyina  the  separation  condition  pf  [q.  f  .  1  i  ana  1* 

••e  are  r, t  non- real  zeros  common  to  all  the  transforms  o(  tne  sea¬ 
ts,  t  n  ei  trie  runst  ion  is  unique  amono  functions  with  the  same  Sup 

t,  n  ow  ever' ,  t»e  funtior.  ’ s  not  necessar  i  I  v  unique  among  ail  f  j< 
ns.  r  i,*-tnerpiore ,  it  was  snown  tnat  if  it  is  known  tnat  one  fare  - 
>.  ac-tter.  f>on  tne  other  by  flipping  at  most  a  finite  number-  c  ‘ 

,  s ,  tti  e”  tne  same  conclusion  nolds  with  the,  condition  cf  Eg.  r .  ’ 

. -3 etc  by  tne  weaker  condition  that  tne  regions  simpl  y  oe  at s  io’r.t 

nenneoat i ve  functions  having  autocorrelations  whose  supports 
’  sty  c^rtam  tr  s  connect  ion  conditions,  it  was  shown  tnat  if  tne 
ri  - f  orms  or  tne  segments  of  the  function  have  no  non-real  zeros  v 
i  ,  t  r  ur  the  fun  ction  is  unique  among  nonneoati  ve  functions. 


ir.aiiy,  it  was  Shown  that  if  the  transform  of  a  function  has 
a  *  i  r,  i  t  e  number  of  non-real  zeros,  tnen  tne  support  of  tne  fur  c 


i  at  well  as  tne  function  itself 

k * i ' "i v*.  led gme n  t  s :  This  WO r  k  was 
■  ier,  t  tic  Research  under  Contra 
d  aener  dl  •  Za  t  i  or  C)f  an  example 


satisfy  some  special  conditions. 

supported  by  the  Air  force  •"'f *  i c 
;t  ho.  F 49620-8C— C-0006 .  "neoren 
communicated  bv  h.  Hurt. 


Append i \ 

p<~Qi|f  of  Tneorem  1:  We  have 

N 

G(w)  =  Vcn  e‘1n6wF(w) 
n=0 


=  F(w) V’cn(e-1Sw)n 
n=0 


N 

=  - (w)aJ'J(e“iSw  -  bn). 

n=l 


e  1  V  h  - 1  new,  ,  , 

H  i  w  1  =  '  3  e  F  ( w  ] 


n=0 


wri  H  r 


n=0 


( w  J  a 


ri  <-*„  -- 


l  Bw 


1) 


nt  d 


H  ie’8"  -  »r 
L  heC 


-  5  f  v,  )c?(  w  ) 


n!  ~  U 

l-r 


/ -E_  e~18w  *  i\ 


-ibw  h 

nEB  \  e  "  bn 


F  or  ur  ft , 


7u 


G  (  u  )  =  G  { u )  . 

’’h  ,j;  , 

H 1  ( u  )  =  F  ( u  )  •  i  G  (  u ) 

=  F(u)  •  G*(u ) 

=  H2(u),  .  Q.E.k 

Proofs  of  Theorems  3  and  4:  The  proofs  of  these  theorems  and 
tr.e  corollary  tc  Theorem  3  use  Corollary  3  to  Theorem  8  and  therefore 
appear  immediately  following  the  proof  of  that  corollary. 
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f  roof  of  Theorem  b:  first,  suppose  rr  =  0.  Ther  ?{(■'•  7 

i’u>  fcL;  (  1\  1,  f  is  intearatle.  Therefore,  by  '3.  ,  p.  ?& 

0  ' 

"heorem  V i , 


FU)  =  F  ( 0 )  e 


- 1  c 


n(’ 


nF(wn) 


hn  } 


1  -  ^ 


nF  wn 


suppose  rs  >  0.  Let 


91  (x)  =  I  f ft)  dt 

Ja, 


=  ,  D  =  bf  aric  c  =  cf ,  Also, 

91  91  T 


r "  .  r 

Q  •,!  >•)  I  ax  =  | 

X 

1  f ( t )  dt  |ax 

'  J  ' 

J 

af  af 

af 

■  r 

X 

1  jf(t )  !  dt  dx 

'af 

'af 

/Df 

|f  (t)  |  dt  dx 

J 

df 

af 

=  (bf  - 

rDf  ■ 

af )  1  Jf (t )  |  dt 

rnerefore  q  is  integrable.  Furthermore 
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'V 


find  let  n  be  tne  inverse  transform  of  H.  To  p-'ov-  .  i  . ,  ’t  . .  <  ♦  v «-  > 
T f  show  that 

;  ■:  n 1  £  : ;  f : . 

_et 


T(w )  = 


w  -  w 

0 


<■  V,  is  sauare-summahle. 


w  w  w 

o  o  o  , 

-  =  —  -  —  !  w  -  v* 

w  -  -  0  0 

—  w  w 

*0  0  o 


if;  v 


H ( w 1  =  -2  c(w'  -  -2  (w_  - 

w  w 

0  0 


w  )  •  '  V. 


-t  ’t  be  tne  inverse  transform  of  T  .  Ther 


t  ' 


lr  /  ,  .  .  i 

*■  sqri f  v  •f  son  i  »  i  j  e 


IVl  > 


wht"'>-  w 


l  .  an  o 

o 


ror  t 


sgn(t)  =  t  0  for  t  =  P 

c  - 1  for  t  <  0 


from  8.3),  we  get 


n(x)  =  ~  f ( x )  -  (wQ  -  wq)(y  *  f )i> 

w  w 

o  o 


To  prove  8.2),  it  suffices  to  show  that 
5.4}  I(Y  *  f)  c  1(f), 
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Nl  >W 


).  > 


,  I  (x-y 

2  I  f  (  V  )[  SQn(  v'r)  -  sqri  (  x  -y}]  e  fy. 


0  ' 


1  r st  suppose  v  ■>  0.  Then 
o 


(Y*f)(xl=--e 


lw  X 

o 


X 

j  f  ( y ) 


-  iw  V 

o-  3 

e  rtv 


'np'-ef  ore,  for  x  <  af,  (y  *  f )( x )  =  0.  for  x 


iw  x  x 


-  iw  y 

<r 


(y  *  f)(x)  =  -  4  e  0  y"f(y)  e  °  dy 


iw  x  r  -iw  y 

;  e  0  y  f (y )  e  0  dy 


IW  X 

c  , 


(w  J 

X' 


=  0. 


>  -t.  y  ’  3  r-  argument  yield?,  tne  same  result  fo  r  <  0.  This  proves 

wr ? c”  implies  8.2)  which  in  turn  implies  8.11. 


13 1  u  -  IX  , u  j 


c  j»- tier  more,  since  O'  u  1  =  cm(ul 


0  X  ue  K. 


G  (  j  )  -  Gm  (  u  )  ‘  =  C(  u  ) c  (  u  )  -  c”  ( u  )  F  (  u  ) 


=  0(  u )  -  c?m  (  U  )  •  r  ( u  ) 

<  4  ,F(u)  l2 

Since  F ( u )  ^  is  summable,  it  follows  from  tne  Lebesoue  dominated 
CO  .mergence  theorem  that 


2 


Rl 


.  o  I  u  )  -  u* 


(u  '  du  - 0 , 

rr>  m-^» 


/ 

'.e.,  o  converges  to  G  in  the  L"  norm,  therefore  converges  to  c 
in  tne  K  norm.  Since 

I q  )  C  I  ( f  i  m  =  1  ,  2 ,  3 ,  .  .  . . 

it  follows  that 

Kg!  c  1(f). 

Now  let.  a(w)  =  6  ( w ) .  We  have 

r,-(w)  =  rip  (w)  -  b(w)  +  e(w). 

'nerefore  0  r  a(w)  <  nr(w)  for  we  (f  &  and  0  z  a(u)  <  2nK  u)  for 
—  —  b  —  ~~  o 

.  Let  2  ■  be  the  distinct  non-zero  zeroes  of  G.  Then 
n  , 
n=l 


a<  Zn ) 


F  (w)  =  G (w) J~[ 


1  -  !L 


n=i  L  "  zn  J 


and  a  similar  argument  yields 


1(f)  £  Kg' 


-oof  of  Theorem  fi:  Assume  ;F(u)  =  G(u)  V-  ue  J\.  First  we 

Her  the  case  N  w  2.  Then  f  * 
the  restriction  of  f  *  f  to  I  _  r 


consider  the  case  N  w  2.  Then  f  *  t  =  g  *  g.  For  n  f  m,  *  f  is 


) f  u  *  q  to  I  -  I  . 
■  n  m 


and 


n  -  and  q  *  q  is  the  restriction 
M  rn  ~n 

Therefore 


f  *  f  =  g  *  gm 
n  m  3n 


F  F  =  G  6  for  n  *  m 
n  m  n  m 
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J' r  i.  t'ir'>v  . ’ ’  S  f  net  it  t^gpr'i  •  '  v 


F  F  F  f-  =  G  G  G  G 
r-  n-  '  ;  '  !  r'?  n?  r;-, 


F  r  =  G  G  . 

ni  n,  ^  n3 

: :  I  '.nest  functions  are  entire,  we  may  divide  E  auction 
;fion  ?.7  and  obtain 


F  F  =  G  G  . 
nn  n„  n„  n, 

c  c  It 


*re‘o re,  ■  *  f.  >2. 


:F  =  G  G  for  n  =  1 . N 

r  n  n  n  ’  ’ 


no’ ds  also  for  N  =  1  since  ff  *  =  GG*  and,  in  t*  <s  case, 
G.  =  G.  Putting  £.5  and  8.8  toqetner; 


:  ■  =  G  G  for  n,  m  s  1,  .  .  !,. 

r .  0  n  If 


jwS  frorr  £.-  tnat 


n< f  i  =  dig  ; 

n  n 


,  r  =  . N.  Now  we  number  tne  elements  nf 


a*  *ne  s^g jt'ire  w 


w,  :  n  =  1  , 
k  ,n 


k  ,  r, 


is  non-decreas inq. 


n  =  l 


’  k  G, 
k  f 


ar 


C  tl  =  ?n.  (0)  =  2:  .  (O’  .  ” her  , 

k  k  k  ’ 


f)  eor  e(T  t, 
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-  1  C  ,  W 

1  10  "l  f 

r  ( w !  =  — F ,  r  t  0  )w  e 
K  Dk  !  k 


r„  'w  ! 
p,  k.  .  r 


n  e  -  e  i 

,  .  k  ,ri  - 
n=  1  '  ' 


k  ,  n  > 


art  0 

-  .  i  ;  1 


, .  ,  .  -1C  w  « 

J  tK  1  ,  ,  Dk  Qk  —  < 
G  i  P  i  w  0 


;  i  f  w  ) 

k  k  .  n 


fl  i  -  , 

,  >  k  ,  n  / 
n=  ■  ’  / 


w  ' 
k  .r  -■ 


Now  from  8.91,  we  obtain 


-  (  W  )  =  ••  ( w  1  +  ( w ' 

rr  "T,  m 


anc;  nence 

■  •  1  "  U  ■  -  •"  !’  W )  =  :  ( W  1  -  n  ( W  1 

r,  ti  tp  m 

*0'-  I  ,  rr.  =  1 ,  .  .  N  and  •¥  wp  (f .  It  follows  from  8.12)  that 
;~*t  not  oe  '  his  equation  is  constant  witn  respect  to  n  and  tne 
i-i  ant  s  •  de  is  constant  with  respect  to  m.  Let 


6  ( w  )  =  n  ( w  1  -  r  ( w ) . 


r,  ft 


e (w i  =  -e(w ) . 


Tnen 


>'  6  (w )  when  g  ( w  )  >  0 

T  {  W  )  =  *> 

[  0  otherwise 


r. 


l  <r  a '  W  }  < 


,  mi  n .  r,  (  w  ) 
l<n<N  n 


from1  wn  i  ch  it  follows  that 
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1  1  -JL 


.  [  Vn  ,  ) 

K  r.  ,  T\ 


,  (  V, 

K  '  K  r ; 


(m  ' 

G  (0! 

k 


(0) 


■  >  "  •  i  0  ,  V  .  i  i  ,  •? .  1 7 ) ,  and  c.lci, 


i  d  w 

G  k  f  v. )  =  A  ^  e  .  ( w ) c  ( v. ' 


1  1  1  VN 


an *:  ^-.20), 


•  (  W-  .  .  1  '  =  1  -V-  vV  E  (£ 


i  w  1  r  f  w  ; 
r  n. 


=  3r,(w)Gm(w* 


lid  -d  )* 

—  nr*  * 

A  A  e-  c  f  w  1 F  ( w ' 

n  ip  n  r 


V 


•ir  . :  V  wf  (f  . 


>v.-'  t.r 


n 


l  r 


wv  or.  t  a 1  n 


ner-  4  **  Qi'  :  .  _  I  ■  an  j  ?  .22  , 


:  0  -9 
?.  fr 


‘  fi  «:■*’*  *  0-‘ 


.  .  . ,  N.  Tneri  t  v  ]  e 


KJV  1& 


.  ' 1  •  ! ■ 

r> 


.  r  , 

=  ‘  V,  ‘  1  /■  ,  ’  = 


C  n  (V;.  t’lflt 


■roir  c.  ,  nt  r.ij  t  *  i  > 


o  v.  .  -  v  C"  v>  .  -  !  v.  . 
!■:  7 ,  it  t.  *  ;  ,ws  that  n  =  c 


iv,  ;  ■.  i"  .  ’  the  ase  r. 


-  Of'!  -  :• 


*  t  - 1 


f  f  '.it:  •  .  i  t  r, 


V 


1  r. 

*  n  ’ 

1  tor 

r  *  ■  . 

* 

n 

* 

*  p 

i'  ‘  r  . 

'  ;■ 

K  t  r '  t-  f '  £  ^  *  1 

r  '  «;  t.  1  ■:  1  r  -i 

f  :  f 

to  :  v.  ■  5 ,  G ,  -  G?o, 

* 

*  * 

♦n-  1 

i r  i 

0*  to  I  , 

.  "np»'ef 

ore,  s 

m:e  ::  =  GG  , 

* 

* 

★ 

c  r  + 

1  • 

F„F  ,  r. 

PM 

r-  r 

O  r;.  U  . 

i  T'  *  0  • "  *-  ,  A  ►“ 

3  I 

v  *  i  j  v  e 

•k 

ik 

*■ 

F  . c  =■ 

u  t  u  -  dn  d 

■  ti 

F  ,  F  „ 

i  i 

=  0,0 
i  2 

and  c.  26',  we  obtain 


15  i IJ  1 


lice  _  r.  ' 

■  '  r  1 


1  -  (6|G,  +  & > F  ^ F  j  •*•  (G ^5, )(&,£, 


1  1  ' 


-  !  F 


F  F  \c  c 

? ?'■  r  l 


★  ★ 

(FlF21tF1F2 


Se'-rt.'re,  '.'net  functions  involved  are  entire,  either 
:  =  G,G*  and,  by  £.25),  r?F*  =  G,3* 


=  G  .C-.-  and,  ay  P.25  I, 


v  c: 


=  g1g1. 


noi as,  tnen  tne  same  argument  used  in  tne  case  N  * 
tTicijsior.  1 ;  of  the  theorem  follows.  I*  ;.,?G'  nolo 


"M  =  G?,  h,  =  6,  and  M  =  H, 


ten 


HR 


m  H  -  F  'f  h  h 

r  ’ n  i  -  r  i  ■ »  r-  ->  ■ 1 : 


★  ★ 


8F  =  FF  . 

■  *  •  •.  *'•  : ■  cj''ih '  t  .j'.eq  ir  the  case  f>  4  ?  applies  ». * t h 
•  ;  '  1  v  a  n,  respectively  and  Conclusion  1  cf 

r  .rows. 

Gov.  assume  to  at  Conclusion  of  tne  tneorem  nolas.  'npn 

’  'Cf-Cn)v' 

6 ( w  )  =  £  c?(w)F  (vw ; 

ft. 

G  ' u '  '  =  G  f  u  :  5  ( u =  c>i  u  )c?  I  u )  f  ( u ) c  !  0  } 

tr 

=  F  (  u  /  F  (  u  ; 

=  F  (  u  I 

■’!  : '  uS  ’ on  "  of  the  theorem  nolds,  then 

1 


"f  '"q 


C  ( w )  F  ( w  ) 


u,  I\, 


G u  )  _  ‘  =  G  (  u  1 G  (o'  =  O  ( u  )c{  u  )  F  f  u  )  F  ( 1. 1 
=  F*(u)F!u' 


f r 00 f  of  Corollary  1  to  Theorem  3:  The  implication 
f  ~  q  Flu},  =  1 G  ( u  1  Vue  J\ 

immed iate . 
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V 


«i»v.  i'  •  ;jMr 


¥  jt  li. 


Since  B  =  t. 


■  ,  if* 

tv.)  =  e 


on  ljv  nv  1'  of  T'iporeir  8  holds,  then 


1  ( c  ,-c 

ie  f  q  , ,  , 

O'  W  !  =  ^  P  -  h  ,  Vv 


q  ( *  i  =  e 1 6  f  (  x  +  cr  -  C  1 
f  c 


lu '•  ior  .  no’  OS ,  then 


-V  ,  19  1!:f+CgW 

u  !  w ;■  =  e  e  y 


F  ;  W  r . 


g  1  * )  =  e 


1  ®  f  /  „  +  -  +  ,  \ 

tf  c  ) 


*  ~  c.  t.r.P 


:-''oc[  >c  '  ro'f  l  1  a '•  v  P  to  "rieoreir  Hece,  N  =  1  and  B 
•:  or  1  1  ar  v  foil  ows .  0. E . P . 


•-•‘00*  ;.»  1 :  By  "•'eonem  5, 


1  ..  <  m  '  ,  -  m 

•  w ,  =  — -  -  ( 0  )  w  e 


■  1  c  ,w 


nF (wn  ’ 


Hi 

n  =  l 


'K^re f c 


r  •  -  r  1  C  rW 

f  ■  1  r  '  ft-  )  I  r,  V  fr‘  T 

(w;  =  —  F  (0)w  e 


n 


nF(wn) 


-v 

n  =  1  n 


rne  lemma  follows.  Q.E.P. 


=  Hf). 


) 
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•  t  oo*  .if  Corol  1  ary  3  to  Tnenrem  Toe  equivalence  of  \  v-. q  "> 
t  immediate. 

')  z>  I*  F ( w )  has  more  than  one  non-real  zero,  then 

fhppinu  one  of  tnem  yields  a  function  G(w)  sucn  that,  by  Coronary 

-  i 

u )  =  f !  u  '  V  ue  I\ 

pi.t  *  and  q  are  not  equivalent. 

1  =£>!;:  By  Lemma  1,  if  f  i  w nas  onl  y  one  non-real  zero  of 

crqer  1,  then  flipping  this  zero  yields  a  function  equivalent  to  f. 

*  is,  in  turn,  equivalent  to  f,  the  implication  follows.  Q.E.P. 

■ ' r oo *  of  Theorem  3:  We  have 

h 

-<  i  -isnwr,  , 

u(w)  =  /  cf)e  F(w) 

n=0 

N 

=  n(w)/  cn(e  ) 

n^O 

N 

=  cfw' (e'18w  -  bn). 
rt=  1 

Let 

M 

O(w)  =  ]~J”  (e“iBW  -  bn) . 

n  =  l 

'ner,  f,  ( w  ;  =  F  ( w  )c7(  w ) .  By  Corollary  3  to  Theorem  8,  since  f  is 
unique,  F  has  either  no  non-real  zeroes  or  one  non-real  zero  of  order 
1.  Therefore,  if  has  no  non-real  zeroes,  then  G  has  either  no  non- 
real  zeroes  or  one  non-real  zero  of  order  1  and,  by  the  same  corol¬ 
lary,  g  is  unique.  Tnerefore,  it  suffices  to  show  that  C  has  no  non- 
1  zeroes.  Suppose  w  )  =  0.  Then  for  some  n  ]  <_  n  _<  N, 
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V 


-  1  t*  u 


- 1  e  Li  e\ 

o  c 


ev 


!ary.  t  ■  ~heorerr  3:  we  have 
r,  r. 


V  ' 


'•w  a^'..  ri 00*' t'rr  j  H'tll 


11 ' 


x  -  e 


iFT, 


r,=  l 


,  M+ 1  . 

:V  =  ?  ,  r,  =  1,  .  .  I,. 

’eerre-T  C  ;  _et  ’ .  =  T  *  g.  Tnen  HI  ft  )  =  ?  '*  }  G f  V.  ). 

<*' :  -  tr  ■'r.-jret  c ,  F  nas  q- tner  no  nor, -real  zeroes  or  me 
z*-<  •  '  ae r  ! .  Sint-.-  r  nas  m  non-real  zeroes ,  h  has  eit 

■'  n i m- r e a  1  z -.•••. jes  or  one  non-real  zero  cf  order  1.  Tnerefore, 
cm  i"  cillery,  i-  ’$  unique.  C.E.D. 

Prm/  of  - r’linia  I' :  The  inequei  it> 

d  <  o 


■  f 


1 s  immediate.  1 t  remains  to  show 


d  j>  u  y  -  a  ^ 
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‘  t  -  a*.  L0t  o  =  T  nen 


g (>*'  * 


fr-' 

J  r 


rry  f *  *■  v  >  ;k  =  ( 


almost  :1  v  >  C.  Let 


at  o 


<->'  ')  =  f  [>'+""a7j ,  -  (x)  =  f(t)f  -  > 


t  -  x  -  a  f »  s-b^.-a^-y 


=  bf  -  af  -  d. 


dll  0 


jc[  t )'.  ( s  -  t 1  dt  =  0 
o 

.a.  s  <•  •  .  i.maer  tnese  conditions,  x  i  tcnmarsh  *  s  fheorem 

p.  ?cs6,  state:'  that  3  *  ana  u  - c?( t )  =  0  a.e.  in  (0,  a}, 
o .  .  t r  !C,  l,  '  and  a  +  u  2  '  ■  Thls  last  ’neauality  implies  tn 
•  :  ?  or  u  >  0  (or  both1,. 

Suppose  *  >  C.  Tner-  fix'.  =  0  a.e.  ir.  (a  s  +  This 
j  the  assumption  tr.at  tne  interval  of  support  of  *  is  *af, 

v  tne  other  hand,  suppose  u  >  0.  Tnen  f(x)  =  0  a.e.  in 
i f  -  u,  b  ) ,  agair  contradicting  tne  interval  of  support  of  * 

iV-Xif  cl  lemna  3:  if  •,  >  d/2,  then  by  Lemma  2 

d  y  >  d  =  b  ^  -  a  ^ 


an  c 


'r.tr.'f  or*. 


<3*  +  Y  ?  b,  -  Y. 


VII  in 
’•  ( * )  =  0 
at 

tonfr  a- 
r.J. 
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\. 


‘t  ’  ^  f  !'f-  =  -  a*,  t- 


■  '•>.  U 

•'*  t»e  i emna  reduce',  the  -  .H_, 


mi*h  assume  >  •  ,•  . 

t’  c  *'•'  ‘->:n3  les t  '  losec  interval  corfa1'  >r>n 

''r  -’it  ’;pst  c  losen  interval  contamvic  r '  * 

-■  :  t;  t’-‘  'lusen  interval  tmitiir.ino  r 

■  -i'-  :  f,  ;x.  ce f  in ec  ty 

1  "  " 

'  ,  =  ,-o,  nt  J 


ui  =  nj. 


J.i  =  ^-d>  -e]. 

■-'c  T  1  a-  and  E  2  d/2.  Furthermore, 

f  *  f  =  0  for  x  c  (  -  6 ,  -a  ly,.,,  f 
1  ■  '  1  '  f  -  a,  define 

i  fix)  for  xe { «f ,  af  +  « ; 

<*  '■ 

I  n  * 

,  0  Other wse 


i;<!  to 

0  Otherwise 


xc.  af  *  a  +  ‘ ,  af 
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Or'  v  .  i  or 


^  (  >■  )q/.  (  >  -  y  )  =  0  V  Xe  1\  . 


•  \?r  vt  •!  ,  b  i, 


/ 


<  /<( 


ri  „  (  x  i  q  t  ;  x  -  y )  dx  <_  /  f  ( x  )  f  (  x  -  y )  dx  =  (f  *  c  (y)  =  0 


rrieref  ore 


r 

J  h{  ( x  )g{  (  x  -  y )  =  0  -V-  yt  K. 


-v  r  i *  i  fj 

*  g{  =  0 

v\»'.  "  cr<  ir.O  ■  1  6:> 

h4  6*  =  0. 

j-.  ( i  it  'or  of  tne  interval  of  support  of  f,  g  t  0.  Therefore 
*  ★  ° 

.  *  arc  G  *  0.  Since  both  H  and  G,  are  entire  functions,  it 

1  6  6 

i  iv.-  tnat  -  =  0  and  h  =  0.  Therefore, 

f  4 

f ! X =  0  for  a. a.  xc^a^  +  a  *  4,  af  +  e]. 

•r.  >  1  was  arbitrary  except  for  tne  condition  0  <  ?  <  B  -  a, 
f ( x '  =  0  f or  a.a .  xe[ a,  +  a,  af  *  b] 

Now  '-edefine  a  and  n  ,  for  0  <  *  <  B  -  a,  by 
6  6 

r  f  (  X  )  for  Xe (bf  -  6 ,  bf ) 

q  ( x  )  = 
y  f  I 

I  0  otherwise 

an  0 


95 


V 


.11'  '  • 


*  ,  _ 

‘  o>~  a .  a  .  «, f  a ,  -  >. ,  r  .  - 

-  t  t 

r  '  — 

■  dv  •  4  -  a,  ", 

--  ’  ■'  ce  G ,  u  ;  =  r  i  ,j !'  y 

p  *  i  _  f  *  ♦ 

r  >  _  t*inr;a 

r  i  '  ••'— 

4.,  a*  4  t  U  fbf  -  |  c 

1  1 

i  c  C  '  a 

r,  ^  -  -  a  u  foc  -  i  c. 

a^.c:  :  .  ;-r-  ■ , 

d,  -  a,  =  o  =  0  -  a 

9  0 

*  =  9 f  v  -  c 

*  a  .  Teen  a.  ^  a*  anc. 

c  n  f  ’ 

3 

ii  i 

1  1 

3.  ,  ^  (i:  u  ( C  f  -  5  c 

tn>>  smallest 

c  1  o S e ."I  interval  contain  inc 

t-ne  smallest 

closec  interval  contain inc 

tri  r-  smallest 

closed  interval  containina 

tne  smallest 

closeG  iriterval  containing 
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u  1  ^ 
u 


.•  *  v  toe  eDa^at'or  cor'd’  tier  r-  .  1  1  .  MirtoeTi. 
••  =  !  ‘  ?  , }  fi  2  !  F  '■  =  t. 

]'.dr\  tc  Theorem  f  ~  h.  Since  also  1  - 

.  o.e.r. 

*Z..±2:  tor  0  *  n'.  *r  =  v  and  r r  ^ 


V, 

w 

r- 


n  1 


l  _ 


w  n 


n ' 


—  I 

K  i  i 

n  '  -J 


‘  n 


L 


p 


J 


i 


ci’  cancellations,  we  may  assume  without  loss 


cur.  to  at  to^  *  a*'e  zeroes  of  t  V 

r. 

*  " u ’ a  :  e>  *.  *'■  show  i no  that 


-  ; :  tc  s 


ar  c  ' 
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G  (  w  =  p  F  (  W  }  i  1  w  )  . 


n=l  n=l 


TT  .Bh 

j  r>  -  v 

n=  1 

oree  o *  the  numerator  of  L  is  less  tnan  tne  dea^ee 

or,  e  i  par  t'i  a  ’  f^dctior  de r  ornp.  '<  '  \  •  Or  9\ 


\  r  n 

V’  N"1  '  n  ,  r 
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we  pause  he-e  to  shov.',  for  later  use,  that 


S.41) 


Cn,6  ^  0  for  n  =  1 . N. 


Let  t  oe  an  integer  1  £  t  <_  N.  We  wish  to  show  that  C^.e^  =  0.  By 
3.^0) , 

-.4-  w  -  wt'  lH(*)  =  *  ctt6t_1(w  -  wt)  *  ct^_2(w  -  wt)- 


+  ct,l(w  -  wt} 


»«■ 


_  1 


N  6n 


(w 


wt)Btrr 

1  (w  -  w. 


rial  r=l  (w  "n ' 

n*t 


By  5.39), 


n<w  -  wn) n  -  ne  -  n 

n=l 


8.43)  (w  -  wt)  tH(w)  =  2=1 - ^ - 

II  (w  ’  wn)  n 

n*l 

h/t 

Equating  the  rignt  sides  of  Equations  8. 42)  and  5.43)  and  setting 
w  =  w  we  obtain 


I>t  -  wn}  n 
n  =  l 


"t ,  8 1  N 


t  0. 


ri(wt  -  wn>  n 

n=l 

n*t 


This  proves  8.41). 
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Nov.  h>  8.39)  anc  S.40) , 


(w) 


ti 


'-II 

n  =  l  r=l 


r 

n,r 


(w  -  wn)r 


i_et 


Tnen  for  real  u,  r  (u)  is  square-summable.  Let  Y  be  its 
■  *  > '  n  ^  f* 

tr ansf orm . 

Tnen 


8.44)  "J  (  X  )  =  ;> - 7 - c-r-r 

n , r  2  (r  -  1 ) ! 


ir  r  1  1WnX 

[sgn  (vn)  -  sgn  {x)]x  e 


wnere  v.  =  u  +  iv  . 
ri  n  n 

We  nave 


anc: 


e„ 

N  n 


«(w)  =  1 


2Z  "n»r  ‘n,rfw)’ 


n=l  r* 1 

G  ( w )  =  p  F  ( w )  1  ( w ) 


N  Sn 


=  PF(W}  -  Cn,rF(w)  rn,r<w)* 


n=l  r=l 


8.45) 


N  Sn 


g  ( x ; 


=  pf(x)  -  PI]Z]  Cn,r(f  *  Vr)(x)- 


n=l  r=l 


inverse 
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-et  T  =  p  n  Q,  .  By  assumption  T  *  (>.  Since  P  and  Q.  are  open 

1,1  ITi  K 

intervals,  T  is  an  open  interval.  For  xcT,  f(x)  =  0  and  g { x )  =  0. 
Therefore,  from  8.45), 


N  6n 


H  L  Cn  ,r ( f  *  =  0  f°r  xeT . 


n=l  r= 1 


(f  *  y  )(x)  = 
v  n ,r  '  ' 


a 

/'<y>Vr( 


x  -  y)  dy 


a 

r[7~~Tyr  _/f(y)[sgn(vn)  +  son ( x  -  y )] 
,r- 1 

,*  -  y)  e  dy. 


irst  assume  vn  >  0.  Then 


.  r  iw  x 

(f  *  Vr)(x)  =  F---TTT  6 


x  f  -i  -  a  vv  y 

n  tt  i/  ,r-l  n-7  „ 

/  f  y  x  -  v  e  dy. 


xtT,  tnen 


f  f  *  y  ) ;  x 1 
n  ,r  ' 


•  r  iw  x  f  ,  -iw  y 

=  -[r  -  iyj  e  dy 

r  iw  x  f  -iw  y 

=  JT~Z  TJ4  e  y)(*  -  y)r  e  dy 


i  1WnX  y  /r  -  l\  s 

~ ryi e  (  s  J  (1x) 


f  r-  ]  c'  'lwn-y 

f_-|  (y)(-iy )  '  e  dy 


r-  1 

iw  x  T — ,  /  , 

.  r  V  l  r  -  1)  r  ('■-1-S 


-1  (ton);ix: 
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For  iicA,  let 


6.49) 


n,s  s 


iS  V1  Cn,r 
IT  (r“T- 


sTT  F-0  1  S)(wn} 


r=s+l 


and  for  ricB,  let 


.50) 


L  tt 


n,s  s 


iS  V1  "h,r  r(r-l 

IT  (r  -  1  -  s)!  h  1 


-1-S) 


(w  ). 

V  n  l 


r=s+l 


»  ,  V 

1W  X 


0L  X  , 


n=l  s=0 

Tnen,  since  a  is  analytic  and  a(x)  =  0  -V-  xeT, 
8.51)  a(x)=0VxeJ\. 

C  laim 


dn  s  =  0  for  0  £  s  <_  6n  -  1  and  n  =  1 . N. 

Proof  of  claim:  Choose  an  integer  t,  1  <  t  <_  N,  -  vt  <_  vn, 
n  s  1,  .  .  .,  N,  ano  if  vn  =  v^.,  then  6n  <_  6^.  Since  a'x)  = 

0  ¥  xc  K, 


n^t 
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=  Tim 


v’ 


£  r^T  °t,< 


_L 

R  s=0 


,s+l 


N  V1 

££v 

rial  S  =  0 
n#t 


1  im 


JL  [ 

,Bt  Jq 


Rs  T(wn-wt)x 
x  e  dx 


N  V 


£  £  an, 

*  L  _  1  _  r\ 


1  im 


n=l  s=0 
n*t 


U 


ft  i  (w 
s  n 
x  e 


rh  js  , 


N  V 


w--e<££ 


n=l  s=0 
n*t 


cL  .  lim 
n*S  0*. 


where 


fP  s  i(Vwt)x 
cn,s(ft)  *  jf  x$e 

i(Vut)R  (Vvn)R 

»  e  e 


TTw-r^y 


£<->»r 


S  ! 


r=l  -  r)  ![1(wn  -  wt)] 


r+T 


r  ft 


s-r 


If  vt  <  vn,  then  it  is  clear  that 


lim  ~T  Dn  s(R)  =  °- 
R  -*»  ®t  n,s 


If  vn  =  V  then  en  <_  et,  s  <_  en  -  i  <_et 


1  and 
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\ 


P 


j_  ps  +  _ s  !ps~r _ 

P  1  n  ll  TTf  (S  -  r)!_Wn  "  wt 


;  w  -  w. 

1  n  t1 


—  0. 


ncefore  P.53,  holds  V-  n  #  t.  It  tnen  follows  from  P.52)  that 
,  =  u.  Now  repeat  the  above  argument  (starting  from  "Proo*  of 


't  ,c  .  -  1 

claim"',  with  repl 
we  ev en tuallv  obtain 


aim"',  with  bt  replaced  by  s  -  1.  By  continuing  this  proceoure. 


.-4  '  a  =  0  for  0  <  s  <  e  -  1  an  a  n  =  1 ,  .  .  . ,  f; . 
h ,  s  —  —  n 

ii’s  completes  proof  of  claim. 

•.en-.  suppose  ncA.  Then  by  8.49)  and  8.54), 

y  77 :  i':-rr;  F)'f1's)<»n>  <-s  <-»„  -  '■ 

r  =  S*  ’ 

. r  r.  mcfn  C  be  defined  by 


n,fn- 


( r,  t  -  7 

’  ii 


P’6n 

h,eri-l 

1 

^  n ,  B 

’  n 

1 

l”n’6n”' 

1 

h,6n-2 

77 

71 

o 


T\  Cn,3  3  !  n,4 


T!  'n,er 
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we  also  define  the  vertical  vector 


V  = 


1 


tl  (-„) 

Then  the  system  of  equations  8.55)  may  be  written  as 

8.56)  C-V  =  0. 

We  wish  to  conclude  that  V  =  0.  It  suffices  to  show  that 

8.57)  det(C)  t  0. 
we  have 


V1 


det<C>  -  L  TT  K,f 


s=0 

By  8.41),  C  =  0  and  therefore  8.57)  holds  and  V  =  0.  That  is, 

n,Bn 


8.58) 


*»  I  ” 

-1  '  n 


wn)  =  0  for  s  =  0 . 6-1 


and  ncA. 

Furthermore,  since  8n  ±  r\p(wn), 

;(S), 


F '  ; ( wn )  --  0  for  s  -  0,  .  .  .,  Bp  -  1 


and  ncA. 

Since  F^s^(wn)  =  F^(wp)  +  F^s^(wn),  we  obtain 
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d.fsyt  F*s1(wn)  =  0  for  s  =  0 . Bp  -  1 

and  tif A. 

From  S . 58 )  and  8.59),  we  conclude  that  for  ncA,  tne  v.  are  zeroes 
oot'i  r_-|  arid  F,  of  orders  >&n. 

If  ncB,  a  similar  argument  yields  tne  same  result.  Therefor- 

all  the  w  ,  n  =  1,  .  .  .,  h,  are  zeroes  of  Doth  F  .  and  F,  of  orde 
n  -  I  l 

—  n ' 

we  Lave 

N 

Gfw)  =  F(w)J~J 
n  =  l 

Define 


"hen  &r(u),  uc  33,  is  square-summable .  Let  gr  be  the  inverse 
transform  of  &r,  r  =  -1,  1.  We  have 


F  =  F_i  Fr 


Therefore , 


arid 

9  =  9.-1  +  9!  . 

how 
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and 


Uf^)  =  [aff 


s 

m 


i 

j 


Kf!) 


=  [t  ,  b,]. 
L  m’  fJ 


By  Theorem  7, 


I(9r)  =  I(fr),  r  «  -1,  i. 


Theref ore , 


Also,  by  (if  necessary)  adjusting  the  values  of  g_-j  and  q,  on  sets 
of  measure  zero. 


9_1 (*) 


g(x )  for  af  <_  x  <_  sm 
I  0  otherwise 


and 


f  g(x)  for  tm  <_  X  <_  bf 


91(x)  = 


0  otherwise  O.E.D. 

Proof  of  Theorem  1 1 :  Let  R ( F )  =  {w:  Im  w  >  0  and  either  weW  or 


wt  w  ■  =  {  w^ :  n  =  1 ,  .  .  . ,  N  ; .  Let 


A  =  D(F)^R(F)  =  {  z  :  n  =  1,  .  .  .,  «}  with  jz  |  <_  z 


n+1 


n  =  1 . 


In  trie  f actor i zati on  for  F  given  in  Theorem  5,  a  finite  number  of  the 
factors  in  the  infinite  product  may  be  reordered  without  affecting 
the  value  of  the  product.  Therefore  we  have 
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no 


and  let  q  be  its  inverse  transform.  Then 


*  f  r  '  ( 0 ) 

r  (w)  =  U21  G(w] 

F  r  (0) 


F(x)  =  ^ttt^  9 ( x  )  a.e. 


F'  ;(0) 


Therefore  S (g )  =  $(f)  and 


S(g ) *  =  s(f )'  =  (-«  ,  -bf)  u  u  (-P  )  u 


(-af,  *) 


Since  0cPm  ,  0e  -  Pm  and 


Therefore,  by  Theorem  10, 


p  n  (-P  )  *  t>. 

rn  m  '  T  v 

o  o 


P  =  -P 
m  m„ 

0  0 


9_1 (x)  = 


g ( x )  for  x  <  0 


9  otherwise 


9i  (x) 


9 { x )  for  x  >  0 


0  otherwise 


(x)  =  -|>ry-0)  g -1  ( x )  a.e. 


Tr)  «i 
F  r  (0)  1 


Al  •  i'. 


r  *  ,  ,  F^^(O  )  f  ,  - 

=  “77T - G^w;. 

F'  ’  ( 0 ) 

8v  Tnft)ceri  10,  eacn  zeW(F)  is  a  zero  of  f^  of  order  _>  srfz;  and 


:  r  o, f  of  "neorerc  12:  oet  R ( F  'i ,  A,  G  ana  g  be  definea  as  i n  tne  r’-Ou 
*  ’t,oererri  11.  Tnen 

a  p  ( z ) 
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Since  Pm  n  (-Pm  )  ?  tf,  it  follows  by  Theorem  10  that.  pm  =  -P^  . 
Therefore,  f  and  a  are  both  zero  a.e.  on  (-t  ,  -s  )  (j  (s^.  t^ ! 


g  i(x)  = 


q(x)  for  x  <  -t 

—  rr 


0  otherwise 


g2(x)  = 


g(x)  for  x  ^  -  V 


0  otherwise 


Let  f  o ( x )  =  f(x)  +  f,(x).  Then 


f ( x )  =  f_i(x)  +  f 2 ( X  ) . 

By  Theorem  10,  each  zeW(F)  is  a  zero  of  F_-j  and  F ^  of  orders  2  °p ( z ) » 


g_ 1 (W)  =  f_1(w)  n 


"1  -  *' 


x  ±  1  -  - 
ZcW(F)  L  z. 


Op(z) 


6o(w)  =  f?(w)  n 


1  -  - 


x  x  h  -  - 

zew(F)  LJ  z. 


Now  let 


g0(*) 


g(x)  for  -s^  <  x  <s^ 


0  otherwise 


T  neref ore 


8.62) 


*  p  r  1  [  0 1 

F  ,(w)  =  G,(w) 

F  r  (0)  ' 


and 


: .  63 ) 


Fo(w)  =  4^  Go(w)- 
0  Fl  '(0) 


Puting  8.60),  8.61),  8.62),  and  8.63)  together,  we  obtain 


°F(z) 


F* ]  (w)  =  F-j  (w)  n 

1  F(r,(0)  1 


'l  -  - 


w 

ZeW(F)  L1  ~  7J 


Q.E.O. 
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APPENDIX  C 

COMMENTS  ON  CLAIMS  CONCERNING  THE  UNIQUENESS 
OE  SOLUTIONS  TO  THE  PHASE  RETRIEVAL  PROBLEM 
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Abstract 

Questions  are  raised  concerning  some  claims  by  A.H.  Greenaway 
and  R.H.T.  Bates  concerning  the  uniqueness  of  solutions  to  the 
phase  retrieval  problem  for  functions  with  disconnected  support. 

A  counterexample  is  presented. 
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pke,ckdIj»o  r*as  blank-not  filmed 


INTRODUCTION 


In  this  letter,  questions  are  raised  concerning  some  claims  by 
Greenaway  [1]  and  Bates  [2].  These  papers  are  concerned  with  the  question 
of  uniqueness  of  solutions  to  the  phase  retrieval  problem.  This  problem, 
in  the  one-dimensional  case,  can  be  stated  as  follows. 

Let  f  be  a  complex-valued  function  on  tne  real  line  which  vanishes 
outside  of  some  finite  interval.  Let  F  be  its  Fourier  transform.  Given 
tne  modulus  of  F  on  the  real  line,  i.e.,  |F(u)|  for  all  real  u,  the  prob¬ 
lem  is  to  reconstruct  the  original  function,  f,  from  this  information. 

The  general  uniqueness  question  is:  Wow  many  otke.fi  ^unctioni,  qff,  zkcsX 
xdick  vani s/i  outstde  05  some  finite.  inX.zn.voJt  and  whose  fcuAcet  traits  ae.ws 
sates 31/  ,G(u)  =  F(u)  1  jjor  ail  azal  u? 

2.  GREENAWAY'S  PAPER  [1] 

Greenaway  considers  a  situation  in  which  the  unknown  function,  f, 
is  known  to  be  zero  outside  of  the  union  of  two  disjoint  intervals  (a,c) 
and  (d,b).  In  other  words 


f  =  g  +  h  . 

where  g  is  zero  outside  of  (a,c)  and  h  is  zero  outside  of  (d,b)  (see 
Fi gure  1 ) . 

Now  let  F,  G  and  H  be  the  Fourier  transforms  of  f,  g  and  h,  re- 
specively,  extended  by  analyticity  into  the  complex  plane: 


00 
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where  w  =  u+iv  and  u  and  v  are  real.  The  modulus  of  F  on  the  real  line, 
i .e .  ,  F(u) i  ,  is  given  . 

The  question  is:  To  what  extent  do  the.  concLiUcm  deicrtbed  above 
dctevncnc  the  -junction  f? 

The  functions 


elaf(x+6)  and  e1af(-x+6), 

where  the  overbar  denotes  complex  conjugation  and  a  and  3  are  real  , 
have  the  same  Fourier  modulus  on  the  real  line  as  does  f.  If  any  of 
these  functions  are  also  zero  outside  of  the  union  of  the  intervals 
!a,c)  and  (d,b),  then  they  satisfy  all  the  requirements  and  qualify  as 
alternate  solutions.  These  solutions  will  be  said  to  be  a* to  clouted 
with  the  solution  f. 

'low  the  revised  question  is:  A/ie  there  ana  other  iotmtioni  not 
xuccuxted  .o-uta  f? 

Let  w  be  a  non-real  zero  of  F,  and  let 

o 


F|  (w) 


F(w) 


w-w 
_ 0 

w-w 

0 


're  function  F^  can  oe  viewed  as  being  gotten  from  F  by  first  remov¬ 
ing  a  zero  at  w^  and  then  adding  a  zero  at  wQ.  In  other  words,  the 

zero  at  w  has  oeen  ’flioped"  about  the  real  line.  Now  for  real  w,w=u, 
o 


i  u-w 

1 _ o 

u-w 

0 
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and  tnerefore 

F^  (u)  =  F(u) 1  for  al 1  real  u  . 

Hofstetter  [3]  and  Walther  [4]  proved  mat  if  f^  is  any  function  which 
vanisnes  outside  of  some  finite  m’erval  ana  F^u)  =  F(u)  for  all  real 
u,  men  F,  is  gotten  from  f  py  fl>  ,ping  arious  sets  of  non-real  zeroes 
of  r  arid  multiplying  py  a  constant  of  modulus  1  and  by  an  exponential 
'•‘unction.  In  particular,  if  F,  is  ootained  from  F  by  flipping  the  set 
of  all  its  non-real  zeroes,  then  its  inverse  transform,  fp  satisfies 

f1 (x)  =  f(-x)  , 

and  tn us  .  if  vanishes  outside  the  union  of  (a,c)  and  (d,b),  then 
f,  is  a  solution  associated  with  f.  (Here,  if  a  zero  of  F  has  multi¬ 
plicity  n,  it  must  be  flipped  n  times.) 

Now  let  2(F)  denote  the  set  of  all  non-real  zeroes  of  F. 

Greenaway  claims  that  if  F^  is  obtained  from  F  by  flipping  any 
■  vivt  subset,  S,  of  Z(F)  (i.e.,  S^Z(F))  and  if  f^  vanishes  outside  of 
tne  union  of  (a,c)  and  (d,b),  then  all  the  points  in  S  are  zeroes  of 
both  G  and  H . 

Thus,  if  G  and  H  have  no  zeroes  in  common  (which  would  usally  be 
the  case  if  g  and  h  are  gotten  more  or  less  randomly  from  the  real 
world),  then  it  would  follow  that  the  only  solutions  are  f  and  its 
ussoci ated  sol  ut ions . 

Greenaway's  claim  is  true  in  the  special  case  in  which  F  has  only 
a  finite  number  of  non-real  zeroes.  (Actually,  Greenaway's  proof  holds 
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only  for  the  more  restricted  case  in  which  F  has  a  finite  number  of  non- 
real  zeroes  of  order  1.  However,  the  case  of  higher  order  zeroes  can  be 
taken  care  of  by  an  extension  of  his  argument.  See  [5].) 

The  following  counterexample  shows  that  Greenaway's  claim  is  not 
true  in  general.  In  this  counterexample,  the  set  Z(F)  is  infinite  and 
S  is  an  infinite  proper  subset  of  Z(F). 

Counterexampl e  : 

Let 


j(x)  = 


1 


l 


1  -  i  x  j 

for 

XI  <  1 

0 

for 

j  X  I  >  1 

See  Figure  2.  Then  the  Fourier  transform,  4>,  of  t  is  given  by 


Note  that  ?  has  no  non-real  zeroes. 
Now  let 


g(x)  =  8c(x) 


and 


h( x }  =  2  ;(x-4)  +  3;(x-8)  +  2;(x-12). 
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Then  3=8L  has  no  non-real  zeroes  and  hence  G  and  H  have  no  non-real 
zeroes  in  common.  Let 

f ( x )  -  g ( x )  +  h (x ) 

=  3; ( x )  +  2;(x-4)  +  3;(x-3)  +  2 :  ( x  - 1 2 ) 


and  let 


a=-l,  c  =  1  ,  d  =  3 ,  b  =  13. 

Lee  Figure  3.  Then  a  c  •  d  •  b,  the  intervals  (a,c)  and  (d,b)  are 
disjoint,  and  f  is  zero  outside  the  union  of  (a,c)  and  (d.b).  The 
Fourier  transform  of  f  is 

F(w)  =  ^3+2e~4lw  i-  3e'8lw  +  2e'12lw]  }(w) 

=  2(e-4iw  +  2)(e’81w  -  ,5e'4’w  ♦  2)  i(w).  (1) 

Now  let 

g-j  (x )  =  4o(x ) 
and 

h^x)  »  7:  (x-4)  +  4:(x-12) . 
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Then  G^=44>  has  no  non-real  zeroes  and  hence  G-j  and  have  no  non-real 
zeroes  in  common.  Let 

f]  (x)  =  g1  (x)  +  h-j  (x ) 

=  4<$(x)  +  7$(x-4)  +  4<j)(x-12). 

See  Figure  3.  Then  is  also  zero  outside  of  the  union  of  (a,c)  and 
(d,b).  The  Fourier  transform  of  f-|  is 


F, (w)  •  (4  ♦  7e"4lw  ♦  4e‘,2l“)  »(w) 

■  2(2e-4i"  *  l)(e'8iw  -  . 5e'4iw  *  2)  S(») 

■  2e-“i»  (e41»  .  2)  (e-8i«  .  ,5e-41"  *  2)  }<»)  (2) 


It  follows  from  (1)  and  (2)  that 


F]  (w)  =  e 


■4iw 


[  e4itV  \ 

\  e"4lw+2  i 


F(w). 
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Now,  for  real  w,w=u, 


L-4iu  e4lu+2  i  , 

7^7  I  =  K 
I  e  +2  i 

There  fore 

i  F-j  (u) ;  F(u)  |  for  all  real  u. 

Thus  f  and  f-j  are  both  solutions  and  it  is  clear  that  they  are  not 
associated. 

In  order  to  see  which  zeroes  must  be  flipped  to  get  F1  from  F,  let 

;  -j  vwj  =  e  +2 

and 

. 2 ( w )  =  e  -  ,5e  +2. 

Then 

F(w)  =  2r1 (w)  T2(w)$(w)  (3) 

and 

F-j  (w)  =  2e'4iw  T^w)  T2(w)*(w)  (4) 
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Since  S  has  no  non-real  zeroes  and  e  is  never  zero,  it  follows  from 
(3)  and  (4)  that 

z(f)  =  z(r1 )  u  z(r2) 

and 

z(F1 )  =  zT?7)  u  z(:2) 

where 

Z(r])  =  •  w:  weZt^)  •  . 

Thus  the  zeroes  of  F  which  are  in  S^Zd^)  are  flipped.  The  sets  Z ( T-, ) 
and  ZC2)  are  given  by 

Z(ri)  = +  \  n  +  \  1og  2:  .±2,. . . | 

and 

Z(r2)  = '  -  Z  tan"1  V^T  +  j  n  +  i  1  log  2 :  n=0  ,+l  ,+2 , . . .  ■  . 

See  Figure  4.  The  flipping  of  the  zeroes  in  S  is  followed  by  multipli¬ 
cation  by  the  exponential  e''ilw.  The  latter  simply  has  the  effect  of 
translating  f^  into  the  proper  position. 
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In  the  above  example  the  function  >  could  be  replaced  by  any  func¬ 
tion  .vnicn  is  zero  outside  of  the  interval  (-1,1)  and  whose  Fourier  trans- 
•orn  no  non-real  zeroes.  For  example,  ;  could  be  replaced  by 

(x)  =  ( p*©)(2x)  , 

1  for  ;  x  J  <1 
0  f  o  r  .  x  '  >_  1 


wnere  *  denotes  convolution,  or  by 


l 


3.  3ATE ' S  PAPER  [2] 

Bates  considers  the  situation  in  which 

N 

f(x)  •  £'„<*> 

n=l 


wnere  each  f  is  zero  outside  of  an  interval  In  and  the  intervals 
I  ,n=l . . .N  are  pairwise  disjoint.  He  claims  that  if  the  Fourier  trans¬ 
forms,  F  ,  have  no  non-real  zeros  common  to  all  of  them,  then  f  and  its 
associated  solutions  are  the  only  functions  with  compact  support  and 
whose  Fourier  transforms  have  the  same  moduli  as  that  of  the  Fourier 
transform  of  f.  Thus  Bates  claims  even  more  than  Greenaway  does. 
Therefore,  the  above  example  is  also  a  counterexample  to  Bates'  claim. 
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A  stronger  separation  condition  on  the  I  which  does  work  C3n  be 
found  in  [5], 
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Figure  Captions 


Figure  1 :  Member  of  the  class  of  functions  with  disconnected 
support.  Note :  Although  the  functions  g  and  h  are 
represented  here  as  positive  real  functions,  they 
can  be  complex-valued. 

Figure  3:  Functions  f(x)  (above)  and  f^(x)  (below)  have  the 
same  Fourier  modulus. 

Figure  4:  Above:  non-real  zerores  of  F.  Below;  non-real 
zeroes  of  F^ .  The  circled  zeroes  are  flipped. 
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APPENDIX  D 


Willi.  IVlrrmininft  the  Support  of  un  Object  from  the  Support 
of  Its  Autocorrelation.*  j  K  kiknup  and  t  k crimmins  Radar 
and  (fpti c  s  Dniston.  Erinrunnu  ntal  Research  Institute  of  Michigan, 
f‘<)  H>> x  *t)l <s.  Ann  Arbor,  Michigan  -1X107.  -In  astronomy,  x-ray 
crystallography  and  other  disciplines,  one  often  wishes  to  reconstruct 
an  object  distribution  from  its  autocorrelation  or,  equivalently,  freon 
the  modulus  of  it^  Fourier  transform  (i.e  ,  the  phase  retrieval  prob¬ 
lem  f  It  is  also  useful  to  be  able  to  reconstruct  just  the  support  of  the 
objec*  tie  .  the  region  on  which  it  is  nonzero).  In  some  cases,  for  ex 
ample,  to  find  the  relative  locations  of  a  number  of  pomtlike  stars,  the 
object's  support  is  the  desired  information.  In  addition,  once  the 
object's  support  is  known,  the  reconstruction  of  the  object  distribution 
by  tb  .  iterative  method1  is  simplified  We  show  several  methods  of 
finding  sets  which  contain  the  support  of  an  object,  based  on  the 
support  of  its  autocorrelation  The  smaller  these  sets  are.  the  more 
information  they  give  about  the  support  of  the  object  Particularly 
small  sets  containing  the  object’s  support  are  given  by  intersections 
of  its  autocorrelation’s  support  with  translates  of  its  autocorrelat  ion’s 
support.  It  will  be  shown  that  for  special  cases  this  gives  rise  to  a 
unique  reconstruction  of  the  support  of  the  object  from  the  support 
of  its  autocorrelation  1 13  min. ) 
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